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Under  steady-state  conditions,  there  are  no  losses  in  superconducting 
wires.  However,  when  subjected  to  alternating  or  transient  magnetic  fields  or 
transport  currents,  losses  in  type-II  superconductors  can  become  significant. 
This  report  deals  with  hysteresis  losses  at  4  K  measured  by  magnetization  and 
complex  magnetic  susceptibility.  The  theoretical  and  experimental  relation¬ 
ships  between  ac  susceptibility  and  magnetization  as  functions  of  dc  field 
were  examined  in  terms  of  the  critical-state  model  as  developed  by  Carr  and 
Clem.  Minor-loop  hysteresis  loss  was  shown  to  be  obtainable  by  direct  measure¬ 
ment  of  loop  area,  from  the  imaginary  component  of  ac  susceptibility,  and  from 
the  reversible  susceptibility  plus  the  critical  current  density  or  full- 
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predictions  of  Minervinl's  two-dimensional  model.  Hysteresis  loss  measurements 
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ratio  of  matrix  material  to  Nb.  A  review  of  internal  fields  in  supercon¬ 
ductors  showed  the  importance  of  demagnetization  factors  and  internal  fields 
for  the  correct  analysis  of  magnetic  data.  A  theoretical  method  of  calibrating 
ac  susceptometers  for  cylindrical  specimens,  which  is  based  on  a  mutual- 
inductance  calculation,  was  developed. 
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TRANSIENT  LOSSES  IN  SUPERCONDUCTORS 
Final  Report 


EXECUTIVE  SUMkARY 

Under  steady-state  conditions,  there  are  no  losses  in  superconducting 
wires.  However,  when  type-II  superconductors  are  subjected  to  alternating  or 
transient  magnetic  fields  or  transport  currents,  losses  can  become  significant. 
Consequently,  they  are  an  important  design  consideration  in  many  practical 
superconductor  applications.  AC  losses  may  be  categorized  into  hysteresis 
losses,  eddy-current  and  coupling  losses,  and  self-field  losses. 

This  report  deals  mostly  with  hysteresis  losses  at  4  K.  The  experimental 
methods  developed  primarily  involve  the  measurement  of  (1)  magnetization  as  a 
function  of  field,  and  (2)  complex  magnetic  susceptibility  as  functions  of  dc 
field  and  ac-field  amplitude  and  frequency.  Most  measurements  were  made  on 
variations  of  commercial  multifilamentary  Nb-Ti  conductors  with  no  transport 
current  in  longitudinal  fields.  However,  some  results  are  reported  for  Nb-Ti 
with  transport  current  and  Nb^Sn  conductors,  both  in  transverse  fields.  In 
general,  the  experimental  relationship  between  complex  ac  susceptibility  and 
magnetization  as  a  function  of  dc  magnetic  field  was  examined  in  terms  of  the 
critical-state  model  as  developed  by  Carr  and  Clem.  Good  agreement  between 
experiment  and  theory  was  found. 

While  the  dc  magnetic  methods  measure  only  the  hysteresis  component  of 
total  loss,  hysteresis  is  usually  the  largest  component  of  loss  below  1  kHz. 
Additionally,  magnetic  hysteresis  measurements  give  information  on  the  remanent 
magnetization  due  to  trapped  flux  in  superconducting  wires.  This  information 
is  useful  in  selecting  conductors  for  specific  applications.  For  these 
reasons,  it  is  likely  that  magnetic  methods  will  emerge  as  a  standard  measure¬ 
ment  technique  for  quantifying  ac  losses. 


In  addition  to  magnetic  techniques,  two  other  ways  to  measure  transient 
losses  in  superconductors  were  evaluated.  These  methods  were  an  electrical 
wattmeter  and  a  calorimeter.  These  two  other  methods  were  found  to  be 
experimentally  rather  difficult  compared  to  the  magnetic  methods.  The 
electrical  method  requires  the  measurement  of  very  small  changes  in  phase 
angle  between  nonsinusoidal  voltage  and  current.  The  calorimetric  method 
requires  the  measurement  of  very  small  incremental  helium  flows  above  the 
normal  static  boil-off.  More  development  would  be  required,  compared  to  the 
magnetic  methods,  for  the  rapid  acquisition  of  useful  materials  information. 


Program  Objectives 


The  following  were  program  objectives  and  their  related  outputs: 


L.  Correlate  hysteretic  losses  in  a  given  material  for  different  sample 
configurations. 


Paper  I:  Hysteretic  losses  were  studied  in  Nb-Ti  superconductors  with 
different  filament  diameters,  numbers  of  filaments,  and  ratios  of  matrix 


material  to  superconductor.  Hysteresis  decreased  with  decreasing  filament 
diameter  and  was  independent  of  conductor  cross-sectional  shape.  The  loss  per 
volume  of  superconductor,  expected  to  be  independent  of  matrix-to- 
superconductor  ratio,  sometimes  showed  significant  unexplained  dependence. 

2.  Study  the  effect  of  matrix  material  on  the  losses  for  a  given  filamentary 
superconductor.  Study  coupling  losses  as  a  function  of  filament  density  and 
filament-matrix  bonding. 

Paper  V:  The  effect  of  local  area  ratio  of  matrix  material  to  superconductor 
was  studied  in  fine-filament  Nb^Sn  conductors.  The  hysteresis  losses  depended 
on  the  interf ilament  separation  and  area  ratio. 

3.  Extend  existing  theoretical  model  for  transient  transport  current  and 
transient  external  field. 

Paper  IV:  The  two-dimensional  model  for  shielding-hysteresis  and  transport- 
current  losses  was  experimentally  verified  for  the  first  time. 

Results 

A  summary  of  significant  results  follows: 

1.  Hysteresis  loss  was  shown  to  be  accurately  predicted  within  a  few  percent, 
using  the  London-Bean  critical-state  model,  from  the  full-penetration  field  as 
a  function  of  maximum  applied  field.  It  was  shown  that  the  full-penetration 
field  could  be  deduced  directly  from  the  hysteresis  loop,  without  reference  to 
the  critical  current  density.  (See  paper  I.) 

2.  Minor-loop  hysteresis  loss  was  shown  to  be  obtainable  in  three  ways:  by 
direct  measurement  of  loop  area,  from  the  imaginary  component  of  ac  suscepti¬ 
bility,  and  from  the  reversible  susceptibility  plus  the  critical  current 
density.  The  loss  expressions  derived  by  Carr  were  shown  to  be  the  same  as 
those  derived  by  Clem  in  the  limit  where  the  reversible  susceptibility  is 
zero.  This  was  the  first  correlation  of  the  different  experimental  magnetic 
methods.  The  relationships  between  complex  magnetic  susceptibility  and 
permeability  were  discussed.  (See  paper  II.) 

3.  A  theoretical  expression  for  the  reversible  susceptibility  plus  experimen¬ 
tal  measurements  of  full-penetration  field  were  found  to  predict  the  measured 
real  and  imaginary  components  of  ac  susceptibility  as  functions  of  applied 
field  using  Clem's  model.  This  was  the  first  test  of  the  applicability  of  the 
model  to  practical  superconductors.  The  full-penetration  field  was  found  to 
be  inversely  proportional  to  applied  longitudinal  field,  thus  following  a 
Kim-type  relationship.  The  shape  of  the  hysteresis  loop  was  used  to  deduce 
the  superconductor  penetration  depth  and  approximate  the  lower  critical  field. 
(See  paper  III.) 

4.  Hysteresis  and  transport  losses  measured  simultaneously  in  a  Nb-Ti  super¬ 
conducting  coil  were  found  to  substantially  agree  with  the  predictions  of 
Minervini's  cylindrical  model.  This  was  the  first  experimental  verification 
of  the  two-dimensional  theory.  (See  paper  IV.) 


5.  Systematic  hysteresis  loss  measurements  In  a  series  of  fine-filament  Nb.Sn 
superconductors  showed  that  the  effective  filament  diameter  is  a  function  of 
interfilament  separation  and  local  area  ratio  of  matrix  material  to  Nb.  The 
effective  filament  diameter  was  obtained  by  applying  the  critical-state  model 
to  measured  hysteresis  loops.  There  would  be  little  filament  Interaction  due 
to  bridging  for  Nb^Sn  interfilament  separation  greater  than  2  pm.  This  was 
the  first  study  of ^filament  coupling  in  Nb^Sn.  The  study  is  expected  to  have 
an  impact  on  the  design  of  practical  wires:  (See  paper  V.) 

6.  A  review  of  internal  fields  in  superconductors  was  prepared.  The  impor¬ 
tance  of  demagnetization  factors  and  internal  fields  for  the  correct  analysis 
of  magnetic  data  was  discussed.  The  review  will  be  a  useful  guide  for  the 
calibration  of  magnetometers  used  in  measuring  hysteresis  losses  in  supercon¬ 
ductors.  (See  paper  VI.) 

7.  A  new  theoretical  method  of  calibrating  ac  susceptometers  for  cylindrical 
specimens  that  is  based  on  a  mutual-inductance  calculation  was  developed.  The 
results  were  in  excellent  agreement  with  calibrations  using  standard  reference 
materials.  (See  paper  VII.) 

8.  A  multifunction  circuit  was  designed  to  protect  instruments  connected  or 
coupled  to  a  superconducting  magnet  against  large  voltage  surges  in  the  event 
of  a  quench.  (See  paper  VIII.) 

9.  A  review  was  prepared  of  the  methods  used  for  characterizing  alloys, 
especially  superconductors,  with  regard  to  matrix-to-superconductor  volume 
ratio,  filament  radius,  density,  and  weight  and  atomic  percents.  This  review 
will  be  useful  for  some  users  and  manufacturers  of  multifilamentary  wires. 

(See  paper  IX.) 

This  project  required  the  development  of  instrumentation  for  magnetic 
measurements  and  computer  programs  for  data  acquisition  and  processing.  It 
resulted  in  the  publication  of  five  papers,  the  preparation  of  four  other  manu¬ 
scripts  (three  of  which  will  be  submitted  for  publication) ,  and  participation 
in  a  workshop  on  problems  in  superconductivity  [Cryogenics  24,  378  (1984)]. 
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Hysteretlc  losses  in  Nb-TI  superconductors 

R.  8.  Goldfarb  and  A.  F.  Clark 

Electromagnetic  Technology  Division,  National  Bureau  of  Standards.  Boulder,  Colorado  80303 

When  subjected  to  transient  magnetic  Helds,  superconductors  exhibit  losses.  At  low  frequencies, 
most  of  the  dissipation  is  hysteretic.  Magnetization  was  measured  in  an  axial  field  for  eight 
multifilamentary  Nb-Ti  superconducting  wires  with  different  filament  sizes  and  different  ratios  of 
copper  to  superconductor.  The  full-penetration  field  Hp  was  estimated  from  the  high-field  ends  of 
the  hysteresis  loops.  The  estimate  of  Hp  provides  a  method  to  assess  the  critical  current  density  Jc . 
There  was  good  agreement  between  measured  losses  and  those  predicted  from  Hp  and  the  peak 
applied  field. 


INTRODUCTION 

An  ideal  type  I  superconductor  is  perfectly  diamagnetic 
below  its  critical  temperature  Tc  and  critical  field  He .  This  is 
a  consequence  of  the  Meissner  effect, 

ff  =  /a*r  +  fWj  =  0.  X  =  M/H=  -  1 ,  (1) 

in  SI  units,  where  B  is  the  magnetic  flux  density,  H  is  the 
magnetic  field  strength,  M  is  the  magnetization,  and  x  is  the 
susceptibility.  The  diamagnetism  arises  from  circulating 
persistent  surface  currents. 

In  type  II  superconductors  the  pure  Meissner  state  ap¬ 
pears  only  below  HcX .  Between  HcX  and  Hcl  (lower  and  up¬ 
per  critical  fields)  the  “mixed”  state  occurs  in  which  surface 
currents  coexist  with  magnetic  vortices,  usually  containing 
one  flux  quantum  each,  formed  by  local  persistent  currents. 
These  vortices  are  analogous  to  paramagnetic  moments.* 1 2 3 4 5 6 7 8 
The  effective  magnetization  (magnetic  moment  per  unit  vol¬ 
ume),  owing  to  surface  and  vortex  currents,  can  be  measured 
in  the  same  way  as  the  magnetization  of  magnetic  materials. 

According  to  the  “critical  state”  model,  as  the  applied 
field  increases  above  Hcl ,  the  magnetic  flux  and  a  current 
density  J  penetrate  toward  the  center  of  the  superconductor. 
The  magnitude  of  J  is  Jc,  the  critical  current  density,  taken 
to  be  independent  of  field.  The  field  at  which  complete  pene¬ 
tration  occurs  is  termed  the  full-penetration  field  Hp .  For  a 
cylindrical  filament  with  a  field  applied  axially, 

H,=Jcr,  (2) 

where  r  is  the  filament  radius. 

In  “hard"  type  II  superconductors,  such  as  Nb-Ti,  pin¬ 
ning  sites  impede  the  motion  of  flux  vortices,  which  gives  a 
higher  value  of  Jc .  Consequently,  when  the  applied  field  is 

TABLE  t.  Physical  properties  of  multifilamentary  Nb-Ti  superconductors. 


cycled,  hysteresis  occurs.  Flux  pinning  is  the  largest  source 
of  bulk  hysteresis,  and  hysteresis  is  the  largest  component  of 
losses  for  ac  fields  below  about  1  kHz.  It  is  important  to 
minimize  these  losses  in  non-dc  superconducting-magnet 
applications  in  fusion,  accelerators,  and  pulsed  power 
sources.  To  minimize  hysteresis,  almost  ail  commercial  Nb- 
Ti  wires  are  composed  of  many  fine  superconductor  fila¬ 
ments,  1  to  25  ft. m  in  radius,  within  a  thermally  stabilizing 
Cu  matrix.  The  weak  diamagnetic  contribution  of  the  Cu  to 
the  effective  magnetization  is  negligible. 

Simple  magnetic  measurements  can  be  useful  to  fabri¬ 
cators  and  users  of  superconductors.  In  this  paper,  measured 
hysteresis  losses  are  compared  to  those  predicted  from  the¬ 
ory  for  a  number  of  superconductors. 


EXPERIMENT 

A  vibrating-sample  magnetometer  (VSM)  was  used  to 
obtain  dc  hysteresis  loops  at  about  4.3  K.  The  magnetic  field 
was  supplied  by  a  superconducting  magnet.  Field  accuracy 
was  estimated  to  be  within  ±  4  kA/m  (50  Oe),  and  magneti¬ 
zation  accuracy  to  be  within  ±  5%. 

Each  sample  consisted  of  about  5-14  1.5-cm-long  seg¬ 
ments  of  insulated  multifilamentary  Nb-Ti  wire,  bundled  to¬ 
gether.  For  such  lengths  the  filaments  may  be  treated  as  un¬ 
coupled.  The  field  and  all  measurements  were  along  the 
longitudinal  axis.  There  was  no  transport  current  applied  to 
the  wires. 

Magnetization  was  calculated  per  volume  of  supercon¬ 
ductor  filaments.  Demagnetization  factor  effects  were  negli¬ 
gible;  the  internal  field  was  taken  to  be  the  same  as  the  exter- 


Sample  Dimensions  without  Number  of  Filament  Cu/Nb-Ti  volume 

number  insulation  |mm)  filaments  radius  I  ami  ratio 

1  d  :=  0.51  180  12.2  146 

2  0.57x0.84  132  13.0  4  98 

3  d  =  0.68  204b  4  3  2.08 

4  d  =  0.64  180  15.0  1.50 

5  0.52x0.70  180  16.4  1.30 

6  0.64  x  0.89  168  199  1.59 

7  d  —  1.01  114  19  5  4  84 

8  0.63  x  0,88  240  19  3  0  91 
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I  dc  hysteresis  loops  for  eight  multifilamentary  Nb-Ti  superconductors.  Note  that  the  vertical  scales  are  not  all  the  same. 


nal  applied  field.  Before  each  set  of  measurements  was  made, 
the  sample  was  warmed  above  Tc  (about  10.2  K  for  Nb-Ti) 
and  then  cooled  in  zero  field. 

MEASUREMENT  OF  HYSTERESIS 

Table  I  summarizes  the  physical  properties  of  the  su¬ 
perconductor  wires  studied.  Figure  1  shows  the  hysteresis 
loops  for  the  samples  for  a  peak  applied  field  H0  of  1.6 
MA/m  (20  kOe).  Occasional  spurious  nonsmooth  behavior 
is  accounted  for  by  slight  deviations  in  temperature  during 
measurement.  The  area  enclosed  by  each  loop  is  equal  to  the 
loss  W, 


W  =  j  H  <j>  M  dH  , 
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where  H  ranges  from  +  H0  to  —  H0. 

The  full-penetration  field  Hp  can  be  estimated  from  the 
hysteresis  loop.  As  the  loop  is  traversed,  the  filaments  go 
from  full  penetration  in  one  direction  to  full  penetration  in 
the  other  direction.  Therefore,  at  a  high-field  end  of  the  loop, 
is  approximately  one-half  the  field  required  to  reverse  the 
magnetization.  This  reversal  field  1HP  is  determined  by  a 
peak  in  d 1 M  /dH 2  and  is  indicated  on  the  loop  for  sample  1 
in  Fig.  1. 

For  the  case  when  the  peak  applied  field  //0  is  greater 
than  Hr ,  Carr  derived  an  expression  for  the  loss,  assuming 
Jc  to  be  independent  of  magnetic  field.2  In  SI  units. 

W  =  [ I  -  Hp/1H„) .  |4| 

Clearly  H0  is  much  larger  than  Hp  for  the  loops  in  Fig.  1. 
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TABLE  II.  Hysteresis  loss  W  for  a  peak  field  H0  of  1.6  MA/m  (20  kOe). 


>, 

V, 


Sample 

number 

Low  from  loop 

area  |kJ/m’| 

(kA/m| 

H, 

(Oe| 

Loss  from 

Hp  |kJ/m}) 

t 

114 

45 

550 

120 

2 

116 

45 

550 

120 

3 

50 

20 

250 

50 

4 

152 

50 

650 

130 

5 

158 

65 

800 

170 

6 

141 

60 

750 

160 

7 

177 

65 

800 

170 

8 

224 

80 

1000 

210 

Losses  calculated  numerically  using  Eq.  (3)  and  computed 
using  Eq.  (4)  are  compared  in  Table  II.  Considering  the  ap¬ 
proximate  nature  of  the  Hp  determination,  the  agreement  is 
good. 

Based  on  Eq.  (2),  Jc  can  be  estimated  from  Hp .  This 
could  be  useful  for  new  superconductor  materials  where  it  is 
difficult  to  measure./,.  directly  in  wire  samples.  In  this  paper, 
Je  is  the  critical  current  density  for  the  filaments,  not  the 
entire  wire.  For  H0>Hp,  Eqs.  (4)  and  (2)  give  Je  in  terms  of 
measured  loop  area  W-. 

Jc  =  iW/HioHor  (5) 

COMPARISON  OF  HYSTERESIS  LOOPS 

Although  different  wires  from  different  manufacturers 
are  listed  in  the  tables,  some  comparisons  may  be  made. 
Equations  (4)  and  (2)  indicate  the  well-known  fact  that  hys¬ 
teresis  decreases  with  decreasing  filament  radius.  This  is  il¬ 
lustrated  in  a  comparison  of  the  losses  in  samples  1  and  4. 
The  wires  are  virtually  identical  except  for  filament  radius. 
An  extreme  case  is  sample  3,  with  a  large  number  of  very  fine 
filaments.  A  comparison  of  samples  4  and  5  shows  that  con¬ 
ductor  geometry  (round  vs.  rectangular  cross  section)  has  no 
significant  effect  on  the  loss. 

Since  magnetization  is  plotted  per  unit  volume  of  super¬ 
conductor  (Nb-Ti)  rather  than  per  unit  wire  volume,  the  cop- 
per-to-superconduc tor  volume  ratio  is  not  expected  to  have  a 
large  effect  on  the  magnetization  values.  A  comparison  of 
samples  1  and  2  shows  little  difference.  However,  samples  6 
and  8,  with  similar  filament  radii,  show  significant  differ¬ 
ences  in  measured  losses. 


The  shapes  of  the  loops  are  somewhat  intriguing.  The 
loops  for  samples  1 ,  2,  3,  and  4  are  generally  trapezoidal  and 
negatively  sloped.  Those  for  samples  5, 6, 7,  and  8  are  gener¬ 
ally  hexagonal  and  more  hysteretic.  The  difference  in  types  is 
highlighted  if  the  initial  branch  of  each  loop  is  compared  to 
the  final  branch  (see  Fig.  1,  sample  I).  For  the  hexagonal  set, 
there  is  an  attempt,  upon  increasing  the  field  from  zero,  to 
retrace  the  peak  of  the  initial  branch.  These  curves  most 
resemble  the  ideal  reversible  magnetization  curve  that  is  ex¬ 
pected  from  the  critical  state  model. 

CONCLUSION 

Measured  hysteresis  losses  for  eight  multifilamentary 
Nb-Ti  superconductors  were  found  to  be  in  good  agreement 
with  predicted  values  based  on  the  maximum  applied  field 
tt 0  and  the  full-penetration  field  Hp .  The  latter,  a  function  of 
critical  current  density  Jc  and  filament  radius  r,  was  estimat¬ 
ed  from  the  hysteresis  loop.  This  estimate  of  Hp  can  be  a 
convenient  method  of  assessing  the  critical  current  density. 
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MAGNETIC  HYSTERESIS  AND  COMPLEX  SUSCEPTIBILITY  AS  MEASURES  OF 
AC  LOSSES  IN  A  MULT I FILAMENTARY  NbTi  SUPERCONDUCTOR 


R.  B.  Goldfarb  and  A.  F.  Clark 
Electromagnetic  Technology  Division 
National  Bureau  of  Standards 
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Magnetization  and  ac  susceptibility  of  a  standard 
NbTi  superconductor  were  measured  as  «.  function  of 
longitudinal  dc  magnetic  field.  The  ac-f ield-amplitude 
and  frequency  dependences  of  the  complex  susceptibility 
are  examined.  The  magnetization  is  related  to  the 
susceptibility  by  means  of  a  theoretical  derivation 
based  on  the  field  dependence  of  the  critical  current 
density.  Hysteresis  losses,  obtained  directly  from  dc 
hysteresis  loops  and  derived  theoretically  from  ac 
susceptibility  and  critical  current  density,  were  In 
reasonable  agreement. 

Introduction 

When  subjected  to  alternating  or  transient  magnetic 
fields,  superconducting  wires  exhibit  losses.  These 
may  be  categorized  into  hysteresis  losses,  eddy-current 
losses,  coupling  losses,  and  self-field  losses.  Bulk 
hysteresis  losses  are  frequency  Independent  and  arise 
mostly  from  flux  pinning  in  the  superconducting  fila¬ 
ments.  There  may  be,  in  addition,  hysteretlc  losses 
due  to  pinning  at  the  filament  surface.  At  high 
frequencies,  eddy-current  losses  occur  within  the 
matrix  and  on  the  wire  surface  (skin  effect).  Coupling 
losses  in  the  matrix.  Important  over  about  1  kHz,  are 
caused  by  the  transfer  of  coupling  currents  between 
filaments.  If  there  is  a  transport  current,  matrix 
eddy-current  losses  arise  from  the  self  field. 

There  are  a  few  general  methods  of  measuring  ac 
losses  in  superconductors.  The  most  direct  i6  calori¬ 
metric,  In  which  a  superconducting  coil  is  exposed  to 
transient  fields  and  the  dissipated  heat  is  measured 
by  the  boil-off  of  liquid  helium.  Another  method  is 
electrical,  by  measuring  the  phase  angle  between  vol¬ 
tage  and  current.  Wattmeters  work  on  this  principle. 

The  final  methods  are  magnetic.  Magnetization 
measurements  are  possible  because  magnetic  moments  are 
associated  with  circulating  persistent  surface  currents 
and,  in  type-11  superconductors  In  the  mixed  state, 
flux  vortices.  There  are  basically  two  types  of 
magnetic  methods.  In  the  first,  magnetlzstlon  (mag¬ 
netic  moment  per  unit  volume)  is  measured  directly 
using,  for  example,  a  vlbratlng-sample  magnetometer. 
Though  lock-in  detection  at  the  frequency  of  vibration 
is  ussd,  this  is  a  dc  measurement,  since  the  time 
required  to  cycle  the  field  is  on  the  order  of  30 
minutes.  Measurements  may  be  made  while  slowly  ramping 
or  stepping  the  field.  The  latter  gives  more  accurate 
results  if  there  Is  any  Imbalance  in  the  pick-up  colls. 

In  the  second  type  of  magnetic  measurement,  a 
sinusoidal  sc  flsld  is  applied  (sometimes  superimposed 
upon  s  dc  field)  end  the  derivative  of  the  magnetiza¬ 
tion  Is  detected  by  induction  in  a  pickup  coll.  If  the 
pickup  coil  is  uncompensated,  the  induced  voltage  Is 
proportional  to  the  permeability  u.  If  the  coll  is 
compensated  (l.e.,  zero  voltage  when  no  sample  Is 
present),  the  Induced  voltage  i9  proportional  to  the 
susceptibility  x •  In  St  units, 
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where  single  primes  indicate  the  real  (dispersive, 
inductive)  components,  the  double  primes  indicate  the 
imaginary  (absorptive,  resistive,  loss)  components,  and 
Uq  Is  the  permeability  of  free  space.  The  Induced 
voltages  may  be  integrated  with  respect  to  field  (or 
with  respect  to  time,  since  field  is  swept  as  a  func¬ 
tion  of  time)  to  yield  the  frequency-dependent 
magnetization. 1 

The  area  enclosed  by  the  dc  or  low-frequency  ac 
hysteresis  loop  of  magnetization  vs.  field  is  a  direct 
measure  of  the  hysteresis  loss.  Alternately,  the 
imaginary  component  of  susceptibility  (or  permeability) 
at  low  frequencies  is  an  indirect  measure  of  the 
hysteresis  loss.  This  paper  will  discuss  dc  hysteresis 
1  oops  and  complex  susceptibility.  The  commonly  used 
method  of  integrated  voltages  will  not  be  discussed. 

The  conductor  carries  no  transport  current. 

Experimental  Methods 

DC  hysteresis  loops  at  about  4.1  K  were  obtained 
using  a  vibrat ing-sample  magnetometer  (VSM) .  Field 
values  up  to  1.6  MA/m  (20  kOe)  were  determined  from  the 
current  in  a  multifilamentary  Nb-Ti  superconducting 
magnet.  No  correction  was  made  for  trapped  flux;  zero- 
field-value  errors  are  thus  estimated  to  be  i4  kA/m 
(50  Oe). 

AC  susceptibility  was  measured  at  4  K  with  a  lock-in 
amplifier  and  a  compensated  susceptometer  of  the  type 
described  in  Ref.  2.  If  the  induced  waveform  is 
nonsinusoidal ,  the  lock-in  amplifier  will  detect  only 
the  fundamental-frequency  component.  Losses  that 
appear  as  harmonics  are  not  detected.  Therefore,  this 
method  is  most  suitable  for  relatively  small  ac  field 
amplitudes.  For  very  small  ac  field  amplitudes,  x"  is 
known  to  be  about  zero,  and  this  fact  is  used  to  adjust 
the  phase  $  on  the  lock-in  amplifier: 

X1  -  Ixfcost,  (2a) 

X”  ■  I  X  I  sin, .  (2b) 

The  phase  setting  was  kept  constant  for  all  fields,  but 
was  readjusted  for  each  frequency.  The  background 
signal  in  the  absence  of  a  sample  was  subtracted  from 
all  measurements.  The  susceptometer  was  positioned  In 
the  bore  of  the  same  superconducting  magnet  as  used 
with  the  VSM.  Both  apparatus  were  computer  controlled. 
Magnetization  and  susceptibility  accuracies  are  esti¬ 
mated  to  be  within  52. 

The  sample  consisted  of  an  insulated,  Nb-Ti  wire, 
cut  Into  14  straight  segments  1 . 5  cm  long,  and  bundled 
together.  At  low  frequencies  and  for  short  samples  (or 
short  twist  pitch),  the  filaments  may  be  taken  to  be 
uncoupled.  The  field  and  all  measurements  were  along 
the  longitudinal  axis.  Before  each  set  of  measurements 
was  made,  the  sample  was  wanned  to  above  its  critical 
temperature  and  then  cooled  in  zero  field.  Demagneti¬ 
zation-factor  corrections  were  negligible;  the  Internal 
field  was  taken  to  be  the  same  as  the  external  applied 
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field.  Magnetization  end  susceptibility  were  calcu- 
lated  per  volume  of  superconductor  filaments.  Mass  and 
volume  measurements  were  made  on  a  length  of  the  wire, 
before  and  after  dissolving  the  Cu  matrix  in  HNO^t  to 
experimentally  determine  the  filament  radius  andJthe 
volume  and  density  of  NbTi.  The  characteristics  of  the 
wire  are  given  in  Table  I.* 

AC  Field  Dependence  of  Complex  Susceptibility 

The  ac  susceptibility  of  the  sample  was  measured  for 
a  few  field  amplitudes  and  frequencies  as  the  dc  bias 
field  was  stepped.  Figures  la-e  show  a  comparison  of 
the  results.  Only  three  branches  are  shown  for  each 
case  to  avoid  redundancy.  The  shapes  of  the  x '  and  x" 
curves  are  dependent  upon  both  field  amplitude  and 
frequency.  The  frequency  dependence  is  small,  which  is 
evidence  that  hysteresis  is  the  major  loss  mechanism. 
For  the  smaller  field  amplitudes  there  are  large 
differences  between  ascending  and  descending  branches, 
and  the  branches  are  not  symetrlcal  for  positive  and 
negative  external  dc  fields  (Figs.  la-c).  These 
effects  may  be  attributed  to  relatively  large  filament- 
surface  losses  when  the  ac  field  is  small.  BusslAre 
and  Clem*  discuss  the  applicable  case  when  losses  are 
associated  with  the  bending  of  vortices  trapped  at  a 
high  angle  to  the  surface,  and  the  field  is  much  less 


than  the  lower  critical  field  H 
facilitate  this  effect. 
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Table  1.  Characteristics  of  NBS  SUM  1457 

Number  of  filaments:  180 
Filament  radius:  12.20  urn 
Copper- to-NbTi  volume  retie:  1.456 
Density  of  NbTi  alloy:  6.004  g/cm* 

»  NbTi-to-lnsulated-wlre  mass  ratio:  0.3025 


External  Field  (kOe) 

-20  -10  0  10  20 


j? 


* 


External  Field  (MA/m) 


External  Field  (MA/m) 


Figure  1.  I»*gln»ry  (x")  snd  real  (x1)  eosiponents  of  ac  susceptibility  at  a  function  of  applied 
longitudinal  dc  field  at  4  K.  AC  field  amplitudes  and  frequenclss  are  shown. 

Erratum:  Page  333,  first  column,  next-to-last  sentence:  "field"  should  be  "ac  field" 
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Figure  2.  Magnetization  M  ve.  dc  field  H  at  4.1  E  obtained  with  a  VSM.  The  Initial 
branch  for  zero  trapped  flux  and  four  nlnor  loopa  (a-d)  are  ehown. 


Theoretical  Derivation  of  Conplex  Susceptibility 

The  following  theoretical  deacrlptlon  of  coaplex 
auaceptlblllty  la  baaed  on  the  davalopaent  by  Claa  for 
the  caae  of  bulk-plnnlng  hyeteratlc  loaaas  In  type- II 
superconductor*. 5  The  aodel  appllaa  whan  the  ac  field 
aaplltuda  la  aaall  and  the  field  aquations  can  be 
linearized.  It  will  be  shown  that  the  dc  hysteresis 
loop  and  the  field  dependence  of  the  critical  current 
density  are  predictive  of  the  ac  susceptibility. 

The  critical  current  as  a  function  of  transverse 
field*  for  the  aaaple  wee  converted  to  critical  current 
density  J  using  ths  total  superconductor  cross- 
sectional  area  (not  Including  aatrlx  aatarlal).  A  good 
fit  was  obtained  to  the  Kla  equation. 


JC(H) 


Jc(0)/(1*H/H0), 


(3) 


for  fields  H  up  to  4.8  MA/a  (60  kOe).  The  adjustable 
paraaeters  J  (0)  and  Hn  were  found  to  be  9.2  CA/a*  and 
1.1  MA/a  (14ck0e).  respectively.  Uhllc  the  longltudl- 
nal-fleld  J  Is  larger  than  the  transverse-field 
J  ,  they  are,  of  course,  Identical  at  zero  field. 
c  The  full-penetration  longitudinal  flald  H  is 
a  function  of  H  owing  to  the  flald  dependence1^!  J  : 


H  (H) 
V 


V- 


where  r  Is  the  fllaaent  radius. 

a  paraaeter  x  as 


x(H)  -  h/H 


h/Jer, 


It  la  useful  to  define 


(5) 


where  h  Is  the  peak  aaplltuda  of  the  ac  flald  that  la 
auperlaposed  upon  H. 

As  shown  by  Cleaand  using  Eqs.  I,  ths  coaplex 
susceptibility  aay  be  expressed  In  terms  of  an 
Idealized  differential  susceptibility  for  type-II 
superconductors  x  s 


where 


where  Mr(y  Is  the  reversible  aagnetlzatlon,  end 


<l+XteV)«i(*)-l. 

(6) 

(1*Xrav)*2(’0, 

(7) 

(8) 

gj(x)  •  x-5xa/16,  0  <  x  <  l,  (9a) 

g ^ (x)  -  [<-l<x-5xa/16)8  ♦  (-4/3x+2-7x/4»13x*/24)slne 

♦  (-l/2*x/2-xa/6)sin26  ♦  (-x/12*xa/24)aln38 

♦  (-x*/l92)sin4a]/e  ♦  1,  x  >  I,  (9b) 


whera  6(x)  *  Isln'^x”1*),  8(1)  •  »,  8(»)  -  0,  and 


g.(x)  •  (4x-2x*)/3», 


0  <  x  <  1, 


g2(x)  -  (4/x-2/x*)/3«,  x  »  1. 


(10s) 

(10b) 


Using  Eqs.  i  and  7,  x'  and  x"  will  ba  calculated  for 
the  aaaple.  The  najor  hysteresis  loop  for  ths  saaple 
la  shown  In  Fig.  2.  The  lower  critical  field  la 
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Figure  3.  Iaaglnary  (x")  and  real  (x')  coapenenta  of 
ac  auaceptlblllty  calculated  for  an  ac  field  aaplltuda 
of  36.2  kA/a  (706  Oe)  using  the  theoretical  aodel,  the 
data  In  Fig.  2,  and  ths  field  dspendanca  of  J£  (Eq.  3). 
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Tabl*  II.  Comparison  of  hyateratic  losses  obtained  by  dc  and  ac  methods 
for  a  field  amplitude  of  56.2  kA/m  (706  Oe) . 


Center  field 


Loss  from  minor- 
loop  area  (kJ/a3) 


Loss  from  x" 
Eg.  11  (kJ/m$) 


Loss  from  Eqs. 
12,  13  (kJ/a3) 


0 

0 

0.6 

0.7 

0.9 

400 

5.0 

1.3 

2.0 

2.1 

800 

10.0 

1.4 

1.8 

2.6 

1540 

19.3 

1.3 

1.6 

2.6 

about  0.12  MA/m  (1.5  kOe) .  The  tangent  to  the  curve 
at  any  H  is  the  differential  susceptibility  xdc.  ** 
an  approximation,  xd  will  be  used  for  Xrcv*  This  is 
somewhat  reasonabledcfor  the  initial  branch  of  the 
hysteresis  loop,  but  rather  unreasonable  for  the  later 
irreversible  branches.  However,  theoretical  x  and 
experimental  xd  are  both  about  zero  for  fields  much 
larger  than  HcjCand  therefore  do  not  greatly  affect  the 
computation. 

A  program  was  written  to  compute,  for  each  value  of 
H  and  xA  »  the  predicted  J  (H) ,  H  (H) ,  x(H) ,  g.(x), 
g2(x) ,  X*(x),  and  xM(*)»  given  r^and  h.  For  h  -  56.2 

kA/m  (706  Oe) ,  x'  and  x”  as  computed  are  plotted  as  a 
function  of  H  in  Fig.  3.  As  in  Fig.  1,  only  the  first 
three  branches  are  shown.  No  adjustable  parameters 
were  used.  A  comparison  with  actual  measured  x'  an<*  XH 
in  Figs.  Id ,e  shows  good  agreement  in  both  shape  and 
magnitude.  Ideally,  the  computed  curves  in  Fig.  3 
would  be  symmetrical  with  respect  to  the  external  dc 
field  and  show  no  difference  between  the  initial 
segment  and  the  return  portion,  as  expected  from 
critical  state  theory.5  Deviations  from  this  are 
because  xdc  was  used  instead  of  Xrev* 

Hysteresis  Loss 

In  this  section,  the  minor  loops  in  Fig.  2  are 
examined.  The  loss  corresponding  to  each  is  calculated 
in  three  ways:  minor-loop  area,  from  x".  and  fro*  Jc- 

The  energy  loss  per  unit  volume  per  cycle  is5 


Substituting  for  x"  using  Eqs.  7,  10b,  and  5  for  full 
field  penetration  (x  >  1), 

W«  (Uxrev)(4u0hJcr/3)(l-Jcr/2h).  (12) 

Except  for  the  factor  (l*x  )»  the  same  equation  Is 
derived,  in  cga  emu,  by  CaiffY7  But  for  H  >>  H^, 

X  a  0.  For  partial  field  penetration  (x  <  l), 
usfng  Eqs.  7,  10a,  and  5,  Eq.  11  may  be  expressed  as 

W-  (^Xrev)(^u0h3/3Jcr)(l-h/2Jcr).  (13) 

again  similar  to  Carr’s  expression  for  this  case. 

Figure  2  shows  four  minor  dc  hysteresis  loops  for 
the  sample  determined  experimentally  with  the  VSM,  The 
peak  amplitude  for  each  minor  loop  Is  56.2  kA/m  (706 
Oe).  The  area  enclosed  by  each  loop  was  calculated 
numerically  by 


W  -  u0  /h  dM  >  u0  /m 


This  hysteresis  loss  Is  compsr.d  in  Table  II  Co  the 
losses  derived  froa  Eq.  11,  using  x"  measured  at  10  Hr 
In  Si. 2  kA/a  (706  Oe)  peak  field  (Fig.  Id).  Another 
comparison  can  be  made  with  losses  obtained  using  Eqs. 
12  and  13,  with  J  (at  ths  center  field)  from  Eq.  3. 
For  the  first  mlno$  loop,  *  <  1.  Once  again,  x 
Is  approximated  by  x,.  from  the  Initial  branch  ofr*v 
che  hysteresis  loop.  The  agreement  Is  reasonsbls. 
Indicating  that  hystsretlc  losses  may  bs  measured 
directly  with  a  hysteresis  loop,  or  estimated  by  mea¬ 
surement  of  ths  Imaginary  component  of  low-frequency  ec 
susceptibility,  or  from  s  reasonable  estimate  of  J  (H) . 


Possible  Inaccuracies  owing  to  ths  following  should 
bs  considered:  (1)  J  (H)  was  measured  In  transverse 
rather  then  longitudinal  field;  (2)  losses  arc  not 
entirely  hystsretlc;  (3)  Induced  voltage  In  suaesp- 
tometer  pickup  coll  contains  higher-order  harmonics; 

(4)  applied  ec  field  amplitude  was  not  smell  enough  to 
allow  linearised  field  equations.  Only  the  first  of 
these  Is  believed  to  be  of  significance,  effecting  the 
data  In  ths  last  column  of  Table  II. 

Conclusion 

Ths  technique  used  to  measure  hystsretlc  losses 
might  depend  on  one’s  application.  If  s  field  is 
stepped  slowly,  a  dc  hysteresis  loop  Is  obtained. 

If  an  alternating  field  Is  swept  and  the  Induced 
voltage  Integrated,  an  ac  hysteresis  loop  may  be 
obtained.  For  the  general  case,  a  dc  bias  field  mey  be 
superimposed  upon  the  stepped  or  alternating  field.  If 
the  Induced  voltage  le  not  Integrated,  but  separated 
into  real  end  Imaginary  parte,  the  complex  suscepti¬ 
bility  is  obtained.  Often,  an  ac  susceptibility 
measurement  le  a  convenient  method.  In  this  pspsr,  it 
was  shown  that  the  complex  susceptibility  In  an  ac 
field  could  be  reasonably  predicted  from  the  dc  hys¬ 
teresis  loop,  the  critical  current  density,  and  the 
filament  radius.  Hysteresis  losses  obtalnsd  directly 
(Eq.  14)  were  compared  to  those  obtained  Indirectly 
using  sxperlmtntal  measurements  of  x"  (Eq.  11)  and 
average  values  of  J£  and  xdc  (Eqs.  12  and  13). 
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ABSTRACT 

DC  magnetization  and  complex  ac  susceptibility  were  measured  at  4  K 
as  functions  of  longitudinal  dc  field  for  a  multifilament ary  Nb-Ti  super¬ 
conductor  with  no  transport  current.  Minor  hysteresis  loops  were  obtained 
in  the  dc  measurements.  The  full-penetration  field,  H  ,  a  function  of 
applied  field,  H,  was  deduced  directly  for  each  minor  lBop.  The  values  for 
H  were  fit  to  the  Kim-type  equation,  H  (H)  •  H  (0) /( 1+H/H^) ,  where  H  (0) 
afid  H.  are  constants.  The  minor  hysterisls-looj)  areas  gave  losses  thSt 
were  in  excellent  agreement  with  Carr's  theoretical  critical-state 
equation,  W  -  (4u-HnH  /3)(1-H  /2Hn) ,  where  Hn  is  the  maximum  applied  field 
for  each  loop.  u  u  p  p  0  0 


An  expression  was  obtained  for  the  ideal  reversible  differential 
susceptibility:  x  "  $q/8ttUq(H-H  ,)A2,  where  is  the  flux  quantum,  H 
is  the  lower  critical  field,  and  AcIs  the  penetration  depth.  H  .  and  A  c 
for  the  sample  were  deduced  from  the  shape  of  the  major  hysteresis  loop. 
Clem's  theoretical  expressions  for  the  real  (x')  and  imaginary  (x") 
components  of  ac  susceptibility  are  functions  of  x  ,  H  ,  and  ac  field 
amplitude,  h.  The  predicted  susceptibilities  based  on  thlse  expressions 
were  in  good  agreement  with  measured  curves  of  x'  and  x"  as  functions  of  h 
and  H.  The  measured  x'  and  x"  were  Independent  of  frequency  up  to  I  kHz, 
as  expected  when  bulk  hysteresis  is  the  primary  loss  mechanism. 


INTRODUCTION 

In  an  earlier  work1  we  discussed  the  relationship  between  dc  magnetiza¬ 
tion  and  ac  susceptibility  in  a  type-II  superconductor.  That  work  described 
how  magnetic  measurements  provide  information  on  hysteresis  losses.  In 
this  paper,  as  in  Ref.  1,  we  examine  magnetization  and  susceptibility  for 
longitudinal  fields  and  no  transport  current.  Here,  however,  we  obtain  the 
full-penetration  fields,  H  ,  directly  from  the  hysteresis  loops  rather  than 
estimate  H  from  measurements  of  critical  current  density,  J  .  Also,  we 
derive  an  Expression  for  the  reversible  susceptibility,  x  ?  rather  than 
using  the  experimental  susceptibility,  xd  >  as  an  approximation.  A  super¬ 
conducting  wire  whose  low-field  magnetization  approached  a  reversible  curve 
was  selected  for  study.2  Several  minor  hysteresis  loops  were  obtained  in 
addition  to  the  major  loop.  The  susceptibility  curves  were  more  nearly 
reversible  and  virtually  independent  of  frequency,  as  expected  from  theory. 
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Table  1.  Characteristics  of  Nb-Tl  Wire 


Cross  section  bare  wire:  0.63  *  0.88  mm 
Twist  length:  1.7  cm 
Cu/Nb-Tl  volume  ratio:  0.91 
Number  of  filaments:  240  » 

Filament  radius:  19.3  um 
Density  of  Nb-Tl  alloy:  6.20  g/cm* 


EXPERIMENT 

The  magnetization  measurements  were  made  with  a  vlbrating-sample 
magnetometer  (VSM)  at  4  X.  Magnetization  was  computed  as  magnetic  moment 
per  unit  volume  of  Nb-Tl.  The  volume  of  the  Cu  matrix  was  not  Included. 

It  Is  useful  to  compare  the  advantages  of  the  VSM  method  vls-b-vls  the 
popular  Integration  method  of  Fietz.*  The  VSM  method  (1)  Is  useful  for 
small  samples.  (2)  Is  a  dc  measurement,  not  sensitive  to  coupling  losses, 
(3)  does  not  require  precise  pick-up  coil  balance  when  the  applied  field  is 
stepped,  and  (4)  Is  not  subject  to  integrator  drift.  The  Integration 
method  (1)  measures  frequency  dependences,  (2)  detects  flux  jumps,  and 
(3)  la  easily  adaptable  to  measurements  with  transport  current. 

The  experimental  methods  are  more  fully  described  in  Ref.  1.  The 
sample  In  this  study  Is  different  and  Its  length  Is  3.0  rather  than  1.5  cm. 
The  purpose  of  a  longer  sample  was  to  avoid  significant  end  effects.  The 
characteristics  of  the  wire  are  given  in  Table  1.  This  is  the  same  wire  as 
sample  8  in  Ref.  2. 


j  DC  HYSTERESIS  LOOP 

Full-Penetration  Field 

The  major  and  minor  hysteresis  loops  are  shown  In  Fig.  1.  The  full- 
penetration  field,  H  ,  a  function  of  applied  field,  H,  may  be  estimated 
I  directly  for  each  mlEor  loop.2  As  the  loops  are  traversed,  the  filaments 

I  go  from  full  penetration  in  one  direction  to  full  penetration  in  the  other 

I  direction.  Therefore,  at  a  high-field  end  of  a  loop,  H  Is  approximately 

one-half  the  field  required  to  reverse  the  magnetizatloE.  The  reversal 
field,  2Hp,  for  the  major  loop  Is  shown  in  Fig.  1. 

The  well-known  Rim  model  for  crltlcsl  current  density  is1* 

Jc(H)  -  ^(Oj/d+H/H^)  ,  (1) 

where  J  (0)  and  H.  are  constants.  In  the  critical-state  model,  for  a  field 
applied  axially. 


where  r  Is  the  filament  radius.  We  did  a  linear  least-squares  fit  of  each 
loop's  Hp  to  the  expression 

Hp(H)  -  Hp(0) / ( 1+H/H^)  ,  (3) 

where  H  (0)  is  a  constant  equal  to  J  (0)r.  The  fit  was  excellent.  We 
obtained "Hp(0)  -  140  kA/m  (1.76  kOe)  Snd  -  1.33  MA/m  (16.7  kOe). 
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Table  2.  Comparison  of  Losses  Obtained 
from  Eq.  4  and  Loop  Area. 


H„ 

H 

Eq. 

4  Area 

(kA/m)  p 

(kJ/ms) 

1591 

63 

165 

159 

1194 

74 

143 

136 

796 

87 

110 

105 

398 

108 

62 

58 

Hysteresis  Loss 

For  the  case  when  the  peak  applied  field,  H-,  Is  greater  than  H  ,  Carr 
derived  an  expression  for  the  loss,  assuming  J  to  be  spatially  conltant 
and  field  Independent  In  the  London-Bean  approximation.5  In  SI  units, 

w*  <^oH0Hp/3)<1-Hp/2H0)  •  (4> 

Clearly  Is  oiuch  larger  than  H  for  the  loops  In  Fig.  1.  Losses  calcu¬ 
lated  by  numerically  Integrating11  the  loops  are  compared  In  Table  2  to 
losses  obtained  from  Eq.  4.  They  agree  to  within  62. 


REVERSIBLE  SUSCEPTIBILITY  x 

rev 

An  exact  equation  for  the  Ideal,  reversible  magnetlzatlon-versus-f leld 
curve  of  a  type-II  superconductor,  H  >  H  ,  Is  not  obtainable  In  analytic 
form.  The  usual  procedure  is  to  divide  &Re  field  Into  three  regions: 
Immediately  above  HcJ  (H  <  H) ,  the  middle  region  (H  <  H  <  H  ) ,  and 
Immediately  below  (Hc<  Hc2) ,  and  make  judicious  approximations. 

One  tractable  case  la  <  H  <  H  _>  when  the  separation  between  flux 
vortices  Is  larger  than  the  coherence  iength,  5,  but  much  smaller  than  the 
penetration  depth,  X.  Thus  the  Glnzburg-Landau  parameter  k  =  X/5  Is  much 
greater  than  1.  For  Nb-Ti  at  4  K,6  X  a  394  nm  and  5*5  nm,  and  this  case 
Is  applicable.  To  first  order,  the  magnetization  as  a  function  of  magnetic 
field  strength  M(H)  is7 

M  -  (♦n/8irynXa)  {ln[4riin<H-H  , ) Xa/*_ ]  +  a}  -  H  .  ,  (5) 


where  $g  is  the  magnetic  flux  quantum  -  h/2e  «  2.068  *  10-15  Wb  (T*ma)  and 
a  is  a  constant  on  the  order  of  unity.  The  functional  dependence  of  M(H) 
can  be  seen  by  dropping  the  a  term  and  using  the  relationship8 


(♦0/4wu0la)  ln(X/5)  . 


M  -  -(♦0/8wW()Xa)  ln[40/4wV()(H-Hcl)5a]  .  (7) 

The  volume  magnetic  susceptibility  x  “  dM/dH  can  be  obtained  from 
either  Eq.  5  or  7  and  is  simply 

X  -  ♦0/8wu0(H-H,1)Xa  .  (8) 

We  take  this  x  to  be  equal  to  the  ideal  reversible  susceptibility  x  in 
the  range  H  ,  <  H  <  H  _.  For  H  <  H  ,  x  -  -1.  X  represents  Be 
susceptibility  under  conditions  of  thermodynamic  equilibrium. 
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Our  goal  la  to  obtain  numerical  values  for  x  for  our  wire.  The 
differential  slope  of  the  third-quadrant  branch  ofe¥he  major  magnetization 
curve  (Fig.  1),  from  0.2  to  1.6  MA/m  (2.5  to  20  kOe),  was  fit  to  Eq.  8 
using  linear  least  squares.  While  this  portion  of  the  magnetization  curve 
Itself  Is  not  reversible,  we  expect  that  its  slope  will  simulate  x  •  For 
the  adjustable  parameters  and  X,  we  obtain  92  kA/m  (1  kOe)  andr^8  nm, 
respectively.  These  values  are  not  too  unreasonable,  though  this  fit  risks 
overestimating  H  ,  with  values  closer  to  H  (0).  Substituting  Into  Eq.  8 
we  get  x  for  ?nis  wire  for  H  .  <  H  <  H  J  (SI  units) : 

T6V  Cl  C4 

Xrev  “  10772/ (H— 9 1 751)  .  (9) 

This  equation  will  be  used  in  the  next  section  to  calculate  x'  end  x". 
An  alternate  approach  would  be  to  use  textbook  values  of  H  .  and  X  in 
Eq.  8. 9  This  would  give  slightly  different  values  of  x'  end  x"  for 
H  t  Hcj.  For  larger  fields,  the  differences  are  negligible. 

MEASURED  AND  IDEAL  x'  AND  x" 

The  real  and  Imaginary  components  of  susceptibility,  x'  end  x"»  were 
measured  as  functions  of  frequency  (10,  100  and  1000  Hz),  sinusoidal  ac 
field  amplitude,  h  [11  and  56  kA/m  (140  and  700  Oe)],  and  dc  bias  field,  H 
[0  -  1.6  MA/m  (0  -  20  kOe) ] .  The  susceptibilities  were  Independent  of  fre¬ 
quency  to  within  3Z.  Since  x"  Is  a  measure  of  the  losses  and  hysteresis  is 
known  to  be  frequency  independent,  bulk  hysteresis  is  probably  the  primary 
loss  mechanism.  Frequency-dependent  eddy-current  and  coupling  losses  are 
likely  insignificant  for  this  sample. 

Figures  2  and  3  show  x'  and  x"  measured  at  10  Hz.  The  initial  and 
decreaslng-fleld  branches  are  plotted.  Some  irreversibility  may  be  seen. 

In  Fig.  3,  h  is  on  the  order  of  H  .  As  shown  by  Clem,10  x'  and  x"  as 

functions  of  H  may  be  predicted  f?om  h,  H  ,  and  x  •  For  our  cases, 

h  <  H  ,  the  equations  are:1,10  p  reV 

P 

X'  “  (l+xrev)(h/Hp-5ha/16Hpa)-l  ,  (10) 

X"  •  (l+Xrev)(4h/Hp-2ha/Hpa)/3ir  .  (11) 

Using  the  expressions  for  H  and  x  (Eqs.  3  and  9),  we  computed  the 
theoretical  curves  of  x'  and  x"  f8r  this  wire  for  ac  field  amplitudes  of  11 
and  56  kA/m  (140  and  700  Oe) .  They  are  shown  in  Figs.  4  and  5.  They  com¬ 
pare  quite  favorably  with  the  actual  curves  in  Figs.  2  and  3.  For  H  <  H  . , 
x'  *>  -1  and  x"  “  0.  There  are  discontinuities  at  H^,  as  expected.  In  £ne 
actual  curves  these  are  rounded  owing  to  flux  pinning. 


CONCLUSION 

DC  magnetization  curves  and  measurements  of  complex  susceptibility 
provide  information  on  hysteretlc  losses  in  multlf ilamentary  superconductors 
such  as  Nb-Tl.  These  measurements  may  be  performed  on  small  samples. 
Interpretation  of  the  data  yields  information  on  full-penetration  field, 
critical  current  density,  lower  critical  field,  and  penetration  depth.  The 
frequency  independence  of  susceptibility  suggests  bulk  hysteresis  as  the 
primary  loss  mechanism. 
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LOSSES  IN  A  Nb-Ti  SUPERCONDUCTOR  AS  FUNCTIONS  OF  AC  FIELD  AMPLITUDE 
AND  DC  TRANSPORT  CURRENT 


M.  Dragomirecky ,  J.V.  Minervini,  J.W,  Ekin,  R.B.  Goldfarb,  and  A.F.  Clark 
National  Bureau  of  Standards,  Boulder,  Colorado  80303,  USA 


Hysteretic  shielding  losses  and  transport  losses  were  measured 
in  a  multifilamentary  Nb-Ti  superconducting  coil  as  functions  of 
transverse  ac  field  amplitude  and  dc  transport  current.  The 
conductor  was  biased  with  a  dc  field.  There  was  significant 
agreement  with  the  predictions  of  Minervini' s  two-dimensional 
theoretical  model. 


INTRODUCTION 

AC  losses  for  the  case  of  current-carrying  conductors  in  transient  applied  magnetic 
fields  are  of  importance  in  practical  applications.  The  theory  for  the  one-dimen¬ 
sional  slab  model  is  well  known  and  a  few  experimental  studies  have  been  made. 
Predominant  among  these  is  the  work  of  Ogaswara,  et  al.  [1-4].  For  multifilamen¬ 
tary  conductors,  they  emphasize  the  dependence  of  loss  on  the  rate  of  change  of  the 
field.  Their  field  amplitudes  were  rather  large,  on  the  order  of  200  kA/m  (2.5 
kOe) ,  Ciazynskl  also  examined  the  dependence  of  hysteresis  and  transport  losses  on 
the  sveep  rate  of  the  applied  field.  The  analysis  was  in  terms  of  an  equivalent 
electrical  network  [5].  Bruzzone  and  Kwasnitza  discussed  the  dependence  of  the 
hysteresis  loss  on  the  ratio  of  ac  field  amplitude  to  fullrpenetration  field  (6]. 
References  [1-6]  are  primarily  concerned  with  coupling  losses  in  multifilamentary 
conductors.  De  Reuver  et  al.  [7]  recently  reported  measurements  of  losses  in 
conductors  carrying  dc  transport  current  in  a  transverse  ac  field.  They  analyzed 
their  data  in  terms  of  the  slab  model  and  found  qualitative  agreement. 

The  first  complete  two-dimensional  analysis  for  a  cylindrical  conductor  with  a 
transport  current  in  a  transverse  ac  field  was  presented  by  Minervini  in  1981  [8]. 
The  theory  is  based  on  the  London-Bean  approximation  where  the  critical  current 
density,  J  ,  is  independent  of  field,  H.  Ignoring  eddy-current  and  coupling  losses 
(for  low-frequency  ac  fields)  ,  the  total  ac  loss  can  be  separated  into  two  com¬ 
ponents.  (1)  Shielding  hysteresis  loss  occurs  in  the  outer  volume  of  the  super¬ 
conductor  filament  not  occupied  by  transport  current.  (2)  Transport  loss  occurs 
owing  to  flux  motion  through  the  central  portion  of  the  superconductor  filament 
when  the  shielding  region  is  saturated  with  shielding  currents.  The  shielding 
hysteresis  loss,  W^,  and  transport  loss,  W{ ,  as  functions  of  transport  current 
density,  J,  are  shown  in  Figs.  1  and  2  [8],  The  total  loss  would  be  equal  to  the 
sum  of  and  W^.  Curves  are  shown  for  different  ac  field  amplitudes,  h.  The 
energies  Wg  and  Wt  are  normalized  to  W_,  J  is  normalized  to  J  ,  and  h  is  normalized 

to  H  (0),  where,  for  filaments  of  diameter  d, 

P 

H  (0)  "  d  J  /  ir  (1) 

p  c 

is  the  full-penetration  field  for  zero  transport  current  (J  *0),  and 

W0  -  U  P0  Hp(0)  '  3  (2) 

is  the  volume  energy  density  when  h  »  H  (0)  and  J  »  0  [9].  H  (J)  is  a  monotoni- 
cally  decreasing  function  of  J.  In  this^paper  we  report  measufements  of  shielding 
hysteresis  loss  and  transport  loss  that  were  made  in  an  effort  to  test  the  two- 
dimensional  model.  The  results  are  in  good  agreement  with  the  model.  More  experi¬ 
mental  data  are  needed  to  adequately  determine  the  detailed  shape  of  the  curves. 
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EXPERIMENT  AND  RESULTS 


A  multifilamentary  Nb-Ti  conductor  was  noninductively  wound  on  a  coil  form.  Inner 
and  outer  balanced  pick-up  coils  were  wound  coaxially.  The  sample  was  immersed  in 
liquid  helium.  It  was  exposed  to  a  3-MA/m  (40-k0e)  dc  bias  field.  The  purpose  of 
the  bias  field  was  primarily  to  operate  on  th^  relatively  flat  portion  of  the  J  - 
versus-H  curve,  where,  as  assumed  in  the  London-Bean  approximation,  J  is  not  a 
strong  function  of  H  [10].  A  secondary  purpose  of  the  bias  field  was  to  reduce  J 
to  a  value  more  accessible  in  the  laboratory  and  to  operate  under  conditions 
similar  to  those  in  practical  applications. 

The  wire  selected  for  study  had  a  diameter  of  0.635  mm,  180  filaments  30  urn  in 
diameter,  and  a  copper-to-superconductor  volume  ratio  of  1.50.  This  wire  was  the 
same  as  sample  4  in  Ref.  [11].  The  sample  was  18  m  long,  wound  in  two  layers  of 
72  and  74  turns.  The  critical  current  while  in  the  bias  field  was  310  A,  giving  a 
filament  J  of  2.44  »  10®  A/m2  at  3  MA/m  (40  kOe) .  H  (0)  was  2.33  *  104  A/m 
(290  Oe)  and  WQ  was  910  J/m3.  p 

Data  were  acquired  with  a  computer-controlled,  four-channel  digital  recording 
oscilloscope.  The  voltage  waveform  of  the  pick-up  coil  was  integrated  with  respect 
to  time  to  give  magnetization,  which  was  plotted  as  a  function  of  the  ac-field 
waveform.  The  resulting  hysteresis  loop  was  numerically  integrated  to  get  the 
shielding  hysteresis  loss.  The  energy  for  the  loss  was  provided  by  the  ac-field 
power  supply. 

The  magnetization  axis  was  calibrated  by  noting  that  the  maximum  magnetization 
value  for  ac  fields  greater  than  H  (0)  for  zero  transport  current  is  2H  (O')  / 3  ]9]. 
An  alternate  calibration  would  be  pto  use  the  initial  susceptibility  inPzero  bias 
field  with  an  ac  field  amplitude  less  than  the  lower  critical  field,  K  ^ .  However, 
the  external  susceptibility  when  the  field  is  applied  along  the  axis  of  the  sample 
coil  depends  on  the  effective  demagnetization  factor,  a  function  of  how  tightly  the 
coil  is  wound  [12],  This  could  not  be  determined  accurately. 

Transport  loss  was  numerically  computed  as  the  integral,  over  two  complete  cycles, 
of  the  product  of  sample  voltage  and  transport  current.  The  energy  for  the  trans¬ 
port  loss  was  provided  by  the  transport-current  power  supply. 

The  frequency  of  the  ac  field  was  1  Hz,  far  below  the  frequency  for  which  coupling 
currents  or  eddy  currents  become  important.  Unlike  in  earlier  studies  [1-7],  the 
frequency  of  the  ac  field  was  not  a  variable.  The  experimental  set-up  is  showr.  in 
Fig.  3.  Particular  care  was  taken  to  avoid  ground  loops  and  thermal  voltages. 

Very  low-noise  (<  3  uV)  preamplifiers  were  used.  Any  dc  offsets  of  the  amplifiers 
were  subtracted  from  the  measured  waveforms.  The  results  are  shown  in  Figs.  4  and  5 

DISCUSSION 

Generally  speaking,  this  was  a  difficult  experiment.  The  agreement  of  the  experi¬ 
mental  results  (Figs.  4  and  5)  with  the  theory  (Figs.  1  and  2)  is  fairly  good.  The 
need  for  more  data  to  fill-in  the  curves  is  noted. 

For  the  case  J  =  0  and  h  >  H  (0),  the  theoretical  expression  for  W  /W  [8] 
reduces  to  p  s 

Wg/WQ  -  2h/Hp(0)  -  1  .  (3) 

Thus,  the  consistency  of  the  values  for  h/H  (0)  and  W  /W  (J  »  0>  in  Fig.  4  can  be 
evaluated  with  reference  to  Fig.  1.  For  exRmple,  thescurve  for  h/H  (0)  -  1  should 
intercept  the  wg/wg  axis  at  unity.  p 

The  theoretical  curves  for  the  one-dimensional  slab  model  are  only  qualitatively 
similar  to  those  for  the  two-dimensional  cylindrical  model.  However,  when  the  data 
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were  normalized  to  the  values  of  H  (0)  and  Wn  for  the  slab  model  and  compared  to 
the  theoretical  slab-model  curves, ^the  quantitative  agreement  was  poor.  For 
shielding  loss,  the  measured  values  were  greater  by  a  factor  of  about  2.  For 
transport  loss,  the  measured  curves  had  a  different  curvature.  Furthermore,  the 
total-loss  curves  (W^  +  W^)  for  our  data  are  in  excellent  quantitative  agreement 
with  the  theoretical  total-loss  curves  for  the  two-dimensional  model  [8]  but  in 
poor  agreement  with  the  theoretical  curves  for  the  slab  model  [1-4,7]. 

Though  our  data  do  not  include  large  ac  field  amplitudes,  we  consider  this  case  for 
the  purpose  of  comparison  to  the  high-field  work  of  others.  When  h/H  (0)  >>  l, 

W  /V  m  2h/H  (0).  In  this  limit,  from  Eqs.  2  and  3,  13 

sup 

Ws  -  8  u0  H  (0)  h  /  3  ,  (4) 

in  agreement  with  Carr's  derivation  [13], 


CONCLUSION 

The  two-dimensional  cylindrical  model  seems  to  be  better  than  the  one-dimensional 
slab  model  for  analyzing  shielding  hysteresis  losses  and  transport  losses  In 
multifilamentarv  conductors  with  round  filaments. 
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Fig.  1  Theoretical  shielding  hysteresis 
loss  for  different  ac  fields  [8]. 
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Fig.  2  Theoretical  transport  loss  for 
different  ac  fields  [8]. 
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HYSTERESIS  LOSSES  IN  FINE-FILAMENT  INTERNAL-TIN  SUPERCONDUCTORS 


R.  B.  Goldfarb  and  J._  W.  Ekin 
Electromagnetic  Technolqgy  Division 
National  Bureau  of  Standards 
Boulder,  Colorado  80303 

ABSTRACT 

Hysteresis  losses  were  measured  on  a  series  of  fine-filament  Nb^Sn 
superconductors  made  by  the  internal-tin  process.  Hysteresis  was  measured  as 
a  function  of  filament  diameter  and  interfilament  separation  using  a  vibrating- 
s ample  magnetometer  in  transverse  field.  Losses  were  greater  than  expected 
from  the  critical-state  model  that  expresses  loss  as  a  function  of  filament 
diameter.  Micrographs  of  the  reacted  wire  cross  sections  showed  some  inter¬ 
filament  bridging  for  all  wires.  This  gave  rise  to  effective  filament  diameters 
that  were  greater  than  actual  diameters.  The  critical  interfilament  separation, 
above  which  the  losses  would  be  expected  to  follow  the  critical-state  model, 
was  determined. 


INTRODUCTION 

Fine-filament  superconductors  are  of  interest  in  ac-loss  studies  because, 
according  to  the  critical-state  model,  hysteresis  loss  is  proportional  to 
filament  diameter.  Thus,  one  would  expect  that  wires  composed  of  fine  fila¬ 
ments  would  exhibit  small  losses.  In  this  study  we  used  a  vibrating-sample 
magnetometer  in  transverse  magnetic  field  to  measure  the  dc  hysteresis  loss 
for  a  series  of  Nb^Sn  superconductors  made  by  the  internal-tin  process.  There 
was  no  transport  current.  Filament  diameter  and  edge-to-edge  interfilament 
separation  were  systematically  varied  to  determine  whether  filament  coupling 
caused  larger-than-expected  hysteresis  losses.  Such  magnified  losses  may  be 


expressed  in  terms  of  an  effective  filament  diameter  that  is  greater  than  the 
actual  filament  diameter.  The  effective  diameter  is  deduced  from  the  hyster¬ 
esis  loss  using  the  critical-state  model. 

There  has  been  some  previous  work  on  hysteresis  losses  in  fine-filament 
superconductors.  Dubots  et  al.  compared  measured  and  calculated  hysteresis  in 
Nb-Ti  conductors  with  filament  diameters  ranging  from  10  to  0.1  ym  [1],  They 
found  that,  for  diameters  smaller  than  1  ym,  effective  filament  diameters 
remained  about  1  ym.  Information  on  interfilament  separation  was  not  given. 
Carr  and  Wagner  attributed  the  asymmetrical  hysteresis  curve  they  observed  in 
Nb-Ti  with  1.6-ym  filament  diameter  to  relatively  large  surface  currents  [2], 
Ghosh  and  Sampson  studied  fine-filament  Nb-Ti  with  filament  diameters  in  the 
range  1-5  ym  [3],  They  note  that,  for  diameters  less  than  3  ym,  the  magnetisa¬ 
tion  (and  therefore  hysteresis  loss)  is  larger  than  expected  from  the  critical- 
state  model.  They  attribute  this  to  degradation  of  the  magnitudes  of  critical 
current  owing  to  filament  necking  or  damage. 

In  this  paper  we  show  that  certain  Nb^Sn  internal-tin  conductors  with 
filament  diameters  smaller  than  5  ym  have  effective  filament  diameters  greater 
than  expected  from  the  critical-state  model.  The  ratio  of  effective  to  actual 
Nb^Sn  filament  diameters  appears  to  be  a  linear  function  of  interfilament 
separation.  The  slope  of  the  line  is  a  function  of  the  local  area  ratio  of 
matrix  material  to  Nb  for  each  wire.  The  greater  effective  filament  diameters 
appear  to  be  a  consequence  of  bridging  between  filaments. 

SAMPLE  PREPARATION  AND  MEASUREMENT 

The  multifilament ary  Nb^Sn  wires  used  in  this  study  were  composed  of 
subelements  containing  150  filaments  each.  The  array  of  subelements  was  sur¬ 
rounded  by  a  tantalum  diffusion  barrier.  Unreacted  Nb^Sn  wires  were  wrapped 


on  the  threads  of  size  10-24  stainless  steel  screws.  The  thread  pitch  was 
such  that  the  windings  were  not  in  contact  with  each  other.  The  screws  had 
been  previously  oxidized  to  prevent  the  wires  from  sticking  to  the  screws. 
The  internal-tin-process  wires  were  reacted  in  vacuum,  with  cold  ends  out  of 
the  furnace,  according  to  the  following  schedule:  1  day  at  340°C,  A  days  at 
580°C,  and  A  days  at  700°C.  This  schedule  is  adequate  to  fully  react  the 
wires  for  the  filament  diameters  in  this  study  [41.  The  residual  tin  content 
of  each  wire  after  reaction  was  estimated  by  the  manufacturer  to  be  8  wt.%. 

After  reaction,  the  wire  coils  were  unscrewed  from  the  stainless  screws 
and  the  long  ends  cut  off.  Each  coil  was  about  1.3  cm  in  height.  The  wire 
volumes  were  measured  using  Archimedes'  principle.  The  wire  diameter  was  the 
same  for  each  sample,  0.681  mm.  Thus,  the  wire  length  for  each  coil  was 
obtained. 

The  average  Nb-filament  diameters  before  reaction  and  number  of  filaments 
in  each  sample  were  provided  by  the  manufacturer.  The  reacted  filament  dia¬ 
meters  and  total  volume  of  Nb^Sn  in  each  coil  were  calculated  based  on  a  pre¬ 
sumed  increase  of  38%  in  filament  cross-sectional  area  after  reaction  [5]. 
Local  area  ratios  of  matrix  to  Nb  were  determined  by  the  manufacturer  from  the 
arrangement  of  the  Nb  rods  in  the  billets  before  drawing:  the  Nb  may  be 
thought  of  as  enclosed  by  adjacent  hexagons  and  the  cross-sectional  areas  of 
Nb  and  matrix  material  measured.  The  average  center-to-center  filament 
distance,  s,  may  be  determined  from  the  local  area  ratio,  L,  and  Nb  filament 
diameter,  d,  by  the  equation 

L  -  (2/3/it)  (s/d)2  -  1  .  (1) 

The  edge-to-edge  Nb  separation  is  simply  s  -  d. 

Table  1  gives  relevant  parameters  for  the  six  samples.  Note  that  Nb 
diameter  plus  Nb  separation  equals  Nb^Sn  diameter  plus  Nb^Sn  separation. 


Both  sums  are  equivalent  to  the  center-to-center  filament  distance  (not 
shown).  The  matrix-to-Nb  local  area  ratios  of  1.5  and  1.8  correspond  to  area 
ratios  for  reacted  Nb^Sn  of  1.0  and  1.2,  respectively. 

Magnetic  moment  was  measured  at  A, 2  t  0.2  K  with  a  vibrating-sample 
magnetometer  in  the  axis  of  the  coil  (i.e.,  field  transverse  to  the  wire  if  we 
ignore  thread  pitch  of  about  5.5°).  Maximum  applied  field  was  1.6  MA/m 
(20  kOe) .  There  was  no  transport  current  and  the  coils  were  open-circuited. 
Magnetization  was  computed  as  magnetic  moment  per  unit  volume  of  Nb^Sn  (i.e., 
not  including  the  volume  of  the  matrix  material) .  Magnetization  was  plotted 
versus  external  (applied)  magnetic  field.  A  typical  magnetic  hysteresis  loop 
for  one  of  the  samples  (sample  18)  is  shown  in  Fig.  1. 

RESULTS 

Hysteresis  loss  was  computed  by,  numerically  integrating  the  loop  area, 

ignoring  the  initial  branch.  According  to  the  critical-state  model,  for 

transverse  field  geometry  [6]: 

W  *  8  Jc  dgff  uQ  HQ  /  3  it  (SI  units)  (2) 

Thus,  an  effective  filament  diameter,  de^,  may  be  defined  that  is  a  function 

of  the  hysteresis  loss,  W,  for  a  maximum  applied  field,  H^.  Since  we  did  not 

have  data  for  the  critical  current  density,  J  ,  at  1.6  MA/m  (20  kOe) ,  we 

c 

deduced  these  J  based  on  the  critical  current  of  each  wire  measured  at  6. A 
c 

MA/m  (80  kOe)  together  with  the  total  Nb^Sn  cross-sectional  area  and  the 

measured  value  of  the  full-penetration  field,  H  ,  at  1.6  MA/m  (20  kOe) .  J  , 

•  P  c 

Hp,  and  filament  diameter,  d,  are  related  according  to  the  critical-state 
model  [6]: 

Hp  ”  d  Jc  /  ir  (SI  units)  (3) 

Hp  may  be  obtained  as  one-half  the  field  required  to  reverse  the  magnetization 


[7].  From  the  hysteresis  loops,  at  1.6  MA/m  (20  kOe)  is  roughly  75  kA/m 
(940  Oe)  for  all  the  loops.  This  suggests  a  factor  of  10  increase  in 
compared  to  measured  Jc  at  6.4  MA/m  (80  kOe).  Applying  this  factor  of  10  to 
each  conductor  gives  the  values  of  Jc  shown  in  Table  2.  We  then  compute  the 
effective  filament  diameters  of  the  wires,  and  additionally,  the  ratio  of  the 
effective  to  actual  Nb^Sn  filament  diameters. 

Optical  and  scanning-electron  microscopy  revealed  that  some  of  the 
filaments  in  each  wire  were  not  entirely  isolated  but  were  connected  by 
"bridges".  An  example  of  this  is  shown  in  Fig.  2.  Such  filaments  had  typi¬ 
cally  lost  their  original  hexagonal  spacing  and  aligned  themselves  in  rows  in 
the  process  of  wire  drawing  [8]. 

DISCUSSION 

As  shown  in  Fig.  3,  the  relationship  between  actual  and  effective  Nb^Sn 
diameters  is  not  monotonic.  However,  some  Insight  is  obtained  by  plotting  the 
ratio  of  effective  to  actual  Nb^Sn  diameters  versus  Nb^Sn  interfilament 
separation  (see  Fig.  4).  It  appears  that  effective  filament  diameter  is  a 
function  of  both  actual  filament  diameter  and  interfilament  separation.  Thus, 
loss  per  volume  of  Nb^Sn  superconductor  is  a  function  of  both  parameters. 

Further  refinement  is  possible  by  distinguishing  among  those  conductors 
with  the  same  matrix-to-Nb  local  area  ratios.  These  are  identified  in  Figs.  3 
and  4  for  the  two  ratios.  We  see  in  Fig.  3  that  there  is  less  of  an  increase 
In  effective  filament  diameter  for  the  conductors  with  the  larger  local  area 
ratio.  However,  we  note  that,  in  Fig.  4,  the  abscissa  (Nb^Sn  interfilament 
separation)  Is  not  Independent  of  local  area  ratio  (see  Eq.  1).  The  apparently 
linear  dependence  of  the  ef fective-to-actual  diameter  ratio  on  interfilament 
separation  has  no  physical  significance. 


1.5  and  1.8,  respectively.  The  critical  separation  is  offered  simply  as  a 
practical  design  parameter  related  to  the  critical-state  model  and  for  guidance 
in  the  manufacture  of  fine-filament  internal-tin  superconductors  such  as 
these. 
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TABLES 


Table  1.  Parameters  for  six  Nb-Sn  multifilamentary  superconductors 


Sample 

Nb 

Diam. 

(pm) 

Nb3Sn 

Diam. 

(pm) 

Local 

Area 

Ratio 

Nb 

Sep. 

(pm) 

Nb-Sn 

Sep. 

(pm) 

No. 

Fil. 

Length 

(cm) 

Vol. 

Nb,Sn 

(mm3) 

3 

4.2 

4.9 

1.5 

2.1 

1.4 

2850 

17.8 

9.71 

6 

4.0 

4.7 

1.8 

2.4 

1.7 

2850 

17.3 

8.55 

9 

3.5 

1.5 

1.5 

1.0 

5550 

17.1 

9.24 

12 

2.8 

3.3 

1.8 

1.7 

1.2 

5550 

8.01 

15 

2.3 

2.7 

1.5 

1.2 

0.8 

9150 

8.89 

18 

2.2 

2.6 

1.8 

1.3 

0.9 

9150 

16.8 

8.04 

Table  2. 


Sample 


Hysteresis  loss  and  effective  filament  dia- 
meters  for  Nb^Sn. 


Estimated  Eff.  Nb„Sn  Eff. /Actual 

J  W  Fil.  Diam.  Diameter 

(GA/m2)  (MJ/m3)  (pm)  Ratio 


Effective  Nb3Sn  Filament  Diameter  (/zm) 


Nb3Sn 


Matrix/Nb 
Area  Ratio 
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A  —  IQ 


Actual  Nb3Sn  Filament  Diameter  (/zm) 


3.  Effective  filament  diameter  as  a  function  of  actual  Nb^Sn  filament 
diameter  for  two  local  area  ratios  of  matrix  to  Nb. 


V 
v/ 


^5? 


V  ••  V  . 

vN>.v« 


.’■r. 


Effective/Actual  Nb3Sn  Diameter  Ratio 


INTERNAL  FIELDS  IN  MAGNETIC  MATERIALS  AND  SUPERCONDUCTORS 
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ABSTRACT 

This  paper  reviews  some  of  the  concepts  needed  for  the  correct  analysis 
of  magnetization  data,  both  for  magnetic  materials  and  superconductors. 
Demagnetization  factors,  initial  susceptibilities,  and  hysteresis  losses  are 
discussed. 


INTRODUCTION 

Magnetization  measurements  have  acquired  an  Important  function  in  the 
measurement  of  ac  losses  in  superconductors.  Many  of  the  traditional  analysis 
methods  used  for  magnetic  materials  apply  to  superconductors  with  the  assump¬ 
tion  of  perfect  diamagnetism.  Some  refinement  is  necessary  for  type-II  and 
multifilamentary  superconductors. 


DEMAGNETIZATION  FACTORS 


The  magnetization  curve  characteristic  of  a  magnetic  material  is  obtained 
by  plotting  magnetization,  M,  versus  internal  field,  H^nt  is  obtained 

by  correcting  the  external  applied  field,  H  hy  the  demagnetization  field, 
^demag  *  w^ere  ®  the  geometric  demagnetization  factor.  In  SI  units. 


H.  =  H  -  DM  . 
int  ext 


(1) 


For  a  sphere,  D  *  1/3.  For  a  long  cylinder,  with  Hgxt  applied  axially,  D  ■  0; 


for  applied  transversely,  D  ■  1/2.  corresponds  to  the  field  H  in 


the  equation  B  =  yrt(H+M).  H.  _  is  sometimes  referred  to  as  the  Maxwell  field. 
U  int 

Here  we  ignore  any  ac  effects  such  as  eddy  currents. 


INITIAL  SUSCEPTIBILITIES 

Often  the  initial  slope  of  the  magnetization  curve  is  used  for 

calibration  of  magnetometers.  We  define  the  internal  susceptibility,  X^nt»  as 

dM/dH.  and  the  external  susceptibility,  y  ,  as  dM/dH  .  Thus,  from  Eq.  1 
int  r  ext  ext  n 

Xint  =  Xext  1  (1  -  °  Xext)  *  (2) 

Xext  =  Xint  1  (l  +  °  Xint)  *  (3) 

Let  us  consider  the  case  of  a  soft  ferromagnet  with  a  very  large  X^nt* 

In  a  small  applied  field  Hext»  we  would  measure  M  such  that,  from  Eq.  3, 

X  „  =  dM/dH  *  1/D. 
ext  ext 

We  can  treat  a  simple  type-I  superconductor  as  a  perfectly  diamagnetic 

material:  X-  ^  =  -1.  In  an  applied  field  H  ,  from  Eq.  3, 
xnt  rr  ext 

W  '  dM/d"ext  '  1/(1>-l)  •  <*> 

Thus,  for  a  sphere,  Xext  =  -3/2.  For  a  long  cylinder  axially,  Xext  =  -1; 

transversely,  Xext  *  -2. 

An  interesting  case  is  that  of  an  infinitely  thin  sheet  of  finite 

dimensions.  For  H  normal  to  the  surface,  D  =  1.  From  Eq.  4,  x  ,  and 
ext  ext 

thus  M,  will  approach  infinity.  While  this  seems  unphysical  at  first,  we  note 
that  M  is  calculated  as  magnetic  moment  per  unit  volume,  V.  M  gets  infinitely 


large  only  as  V  approaches  zero.  Furthermore,  for  >  0,  flux  immediately 

begins  to  penetrate  the  sheet  in  this  "intermediate"  state;  full  penetration 


occurs  at  H  *  H  ,  the  critical  field, 
c 

For  these  cases,  we  plot  M  versus  H  t  in  Fig.  1  [1],  Note  that  M  versus 

H^n  would  be  characteristic  of  the  material,  independent  of  geometry.  The 

plot  of  M  versus  H,  would  be  the  same  as  that  for  the  axial  cylinder  (D  =  0) 
int 
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The  same  initial  X^nt  would  be  obtained  for  a  type-II  superconductor. 

The  lower  critical  field,  H  , ,  is  signaled  when  the  initial  M-versus-H. 

c  I  int 

curve  deviates  from  linearity,  indicating  flux  penetration  and  the  "mixed" 
state.  (The  intermediate  state  is  also  possible  in  type-II  superconductors.) 


HYSTERESIS  LOSS 

We  now  consider  hysteresis  loss  in  superconductors.  Theoretical  loss 
expressions  are  usually  given  in  terms  of  the  applied  field,  H  .  While  a 
hysteresis  loop  plotted  as  M  versus  is  generally  shaped  differently  from 

one  plotted  as  M  versus  H  the  areas  enclosed  by  the  loops  are  the  same. 
This  can  be  seen  by  integrating  as  a  function  of  M  and  referring  to  Eq. 

1,  as  we  now  show. 

Consider  a  hysteresis  loop  plotted  as  M  versus  H  .  Let  the  end  points 
of  the  loop  be  denoted  as  1  and  2.  The  two  points  are  joined  by  curves  A  and 
B.  Using  the  usual  convention  of  integrating  H  as  a  function  of  M,  the  area 
of  the  loop  enclosed  by  curves  A  and  B  is 


2 

W  =  /i 

H? 

int 

dM  - 

2 

f\ 

H® 

int 

dM 

(5) 

2 

=  / 1 

ha 

ext 

dM  - 

D 

2 

/ 1  M 

dM  - 

2  R  2 

fl  dM  +  D  / 1  M  dM 

1  ext  1 

2 

=  h 

ha 

ext 

dM  - 

2 

/ 1 

hb 

ext 

dM  . 

(6) 

Thus,  the  loop  areas  are  the  same  for  both  M  versus  H.  t  and  M  versus  H 

v  xnt  ext 


MULTIFILAMENTARY  SUPERCONDUCTORS 

The  fields  in  the  somewhat  special  case  of  a  multifilamentary  super¬ 
conductor  have  been  analyzed  by  Carr  [2-4].  In  the  equations,  M  now  refers  to 
the  magnetic  moment  per  unit  volume  of  superconductor,  not  including  matrix 
material.  Often,  however,  M  is  expressed  as  magnetic  moment  per  total 


volume  and  the  volume  fraction  of  superconductor  is  carried  in  the  equations. 
For  hysteresis  calculations,  the  local  field,  applied  to  each  filament 

is  used  rather  than  applied  to  the  composite  conductor.  H^oc  Is  some¬ 

times  referred  to  as  the  Lorentz  or  cavity  field: 

H.  =  H  +  DM  .  (7) 

loc  int 

For  H  applied  axially  to  the  conductor  (and  filaments),  H,  =  H  =  H, 

ext  3  loc  ext  int 

For  Hext  applied  transversely,  +  M/2.  For  the  special 

case  when  is  changed  slowly  and  there  is  no  transport  current,  Eq.  1 

applies:  H,  =  H  -  M/2.  Thus,  Hn  =  H  .  Therefore,  we  can  continue  to 

rr  int  ext  loc  ext 

use  in  the  hysteresis  loss  calculations  for  this  special  case.  Similarly, 

X^oc»  equal  to  dM/dH^^  for  multifilamentary  conductors,  will  equal  x  For 

the  general  transverse  case,  however,  shielding  effects  prevent  us  from  using 
H^c  and  interchangeably. 

We  summarize  the  equations  for  Xint  and  Xext  when  there  is  no  shielding 
and  the  volume  fraction  of  superconductor,  X,  is  carried  in  the  equations 
[2,4].  We  have,  for  applied  axially, 

xi»t  ■  w  ■  -1  •  (8) 

The  expressions  for  Hex<_  applied  transversely  are 

Xlnt  *  -2X/(1+X)  ,  (9) 

W  '  '  <10) 

Equations  9  and  10  are  related  by  Eqs.  2  and  3. 


EFFECTIVE  DEMAGNETIZATION  FACTOR  FOR  COILS 
Sumiyoshi  et  al.  derived  an  effective  demagnetization  factor,  for  a 

single-layer  coil  of  multifilamentary  superconductor  parallel  to  the  coil  axis 
[5].  That  is,  Hgxt  is  applied  transversely  to  the  conductor. 

D  «  [1  -  cos(irr/d)]  /  [cosh(irr/d)  -  cos(irr/d)]  (11) 
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Here,  r  is  the  radius  of  the  filament  bundle  and  2d  is  the  center-to-center 
separation  of  the  bundles.  ^eff  takes  into  account  the  interaction  between 
the  turns  of  the  coil.  We  assume  that  there-  is  no  current  in  the  conductor. 

In  the  limiting  case  where  d  >>  r,  De^  **  1/2,  as  for  Isolated  transverse 
cylinders.  In  the  other  limit,  d  *■  r,  De^  *  0.1588.  Thus  from  Eq.  4,  Xext 
can  range  from  -2  to  -1.189. 

The  effect  of  coil-turn  interaction  has  also  been  discussed  by  Carr  et 
al,  [6],  Zenkevitch  and  Romanyuk  [7],  and  Campbell  [8]. 

CONCLUSION 

The  correct  analysis  of  magnetization  data  requires  that  one  account  for 
demagnetization  factors  and  Internal  fields.  This  applies  to  magnetic  materials 
as  well  as  superconductors. 
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Figure  1.  Magnetization  versus  external  applied  field  for  an  ideal  type-I 

superconductor  of  various  geometries.  The  external  susceptibili¬ 
ties  are  given  by  the  initial  slopes.  The  broken  lines  indicate 
flux  penetration  for  the  intermediate  states.  After  Lynton  [1]. 
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Calibration  of  ac  suacaptometer  for  cylindrical  spacimena 

R.  8.  Goldfarb  and  J.  V.  Minervini 

Electromagnetic  Technology  Division,  National  Bureau  of  Standards,  Boulder,  Colorado  80)03 
(Received  1 1  November  1983;  accepted  for  publication  30  January  1984) 

The  absolute  magnetic  susceptibility  of  cylindrical  specimens  is  obtained  with  an  ac 
susceptometer  whose  calibration  is  based  on  a  calculation  of  mutual  inductance.  An  axially 
magnetized  cylinder  is  modeled  as  a  solenoid  of  the  same  size.  The  mutual  inductance  between 
such  a  solenoid  and  a  pickup  coil  of  arbitrary  dimensions  is  computed.  The  susceptibility  is  then  a 
function  of  the  mutual  inductance,  the  cylinder  length,  the  magnitude  and  frequency  of  the  ac 
magnetizing  field,  and  the  voltage  induced  on  the  pickup  coil.  Demagnetization  factor  and  eddy- 
current  effects  are  considered,  an  example  is  given,  and  pickup  coil  compensation  is  discussed. 
Other  calibration  methods  are  also  presented. 

PACS  numbers:  07.55.  +  x 


INTRODUCTION 

A  useful  method  of  measuring  the  magnetic  susceptibility  of 
nonferromagnetic  materials  employs  an  ac  field  source  and  a 
pickup  coil.  A  typical  circuit  is  shown  in  Fig.  1 .  Such  pickup 
coils  usually  require  calibration  with  standards  of  similar 
size  to  that  of  the  specimens  under  study,  except  for  the  case 
of  thin  pickup  coils  with  ellipsoidal  samples,  whose  dipolar 
fields  allow  a  closed-form  mathematical  solution. 1 

A  method  of  calibration  of  an  ac  susceptometer  for  rela¬ 
tively  long  cylindrical  specimens  of  arbitrary  dimensions, 
coaxial  with  a  pickup  coil,  is  based  on  a  calculation  of  mutual 
inductance.  The  mutual  inductance  between  two  coaxial 
coils  (inner  and  outer)  may  be  numerically  computed.  Since 
the  flux  density  distribution  of  a  uniformly,  longitudinally 
magnetized  cylinder  of  uniform  composition  is  the  same  as 
that  of  a  solenoid  of  the  same  size,  a  cylindrical  specimen 
may  be  modeled  as  an  inner  coil.  The  calculation  of  mutual 
inductance  leads  to  the  absolute  magnetic  volume  suscepti¬ 
bility  of  the  cylinder,  when  an  ac  field  is  applied,  in  terms  of 
the  emf  induced  on  the  outer  pickup  coil. 

I.  CALIBRATION  BY  CALCULATION  OF  MUTUAL 
INDUCTANCE 

A.  Equivalent  solanoM 

In  general,  the  flux  density  B  may  be  thought  of  as  due 
to  Amperian  currents  of  volume  density  V  x  M  and  surface 
density  Mxfi,  where  M  is  the  magnetization  and  h  is  the 
unit  vector  normal  to  the  surface.  The  flux  density  B  external 
to  a  cylinder  of  uniform  axial  magnetization  (VxM  =  0)  is 
equivalent  to  that  of  a  model  solenoid  of  the  same  size  with 
ampere-turn  density  NI  /l,  where  M  has  units  A/m,  N  is  the 
number  of  turns,  I  is  the  current,  and  /  is  the  length: 

M  =  NI  /I.  (1) 

The  near  field  cannot  be  expressed  in  closed  form,  but  may 
be  computed  by  superposition  of  B  due  to  elemental  uni¬ 
formly  magnetized  disks  or  elemental  circular  current  loops 
for  the  cylinder  or  equivalent  solenoid,  respectively.  The 
near  field  of  a  current  loop  is  well  known.2-3 

7*1 


The  total  magnetic  flux  (flux  linkage)  0  sensed  by  the 
pickup  coil  may  be  calculated  from  either  B  or  the  vector 
potential  A: 

0  =  B-da  =  <j)A-ds,  (2) 

where  da  is  the  incremental  area,  ds  is  the  incremental  con¬ 
tour,  and  the  integrations  are  over  the  entire  pickup  coil.  The 
mutual  inductance  per  solenoid  turn  is 

L*  =  L,2/N=0/NI=  -i-^A-ds,  (3) 

Where  L ,  2  is  the  mutual  inductance  between  the  solenoid  and 
the  pickup  coil.  L ,  2  is  proportional  to  N\  L  •  is  independent  of 
N  but  depends  on  the  geometries  of  the  solenoid  (i.e.,  speci¬ 
men)  and  the  pickup  coil,  and  the  number  of  turns  on  the 
pickup  coil.  A  general  mutual  inductance  calculation  is  de¬ 
tailed  in  Sec.  I  B.  The  emf  induced  in  the  pickup  coil  due  to 
the  solenoid  is 

v  =  -  d0 /dl  =  ~LndI/dt.  (4) 


Fig.  I .  One  method  [see  number  (6)  in  Tebie  II]  of  measuring  ac  susceptibil¬ 
ity. 
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Converting  current  to  magnetization  using  Eq.  1, 

v  =  -L*ldM/dt.  (5) 

To  this  point,  the  ac  field  source  has  not  been  part  of  the 
derivation.  In  fact,  it  has  been  and  will  continue  to  be  as¬ 
sumed  that  the  emf  induced  in  the  pickup  coil  by  the  field 
source  is  nulled  (see  Sec.  II  B).  The  ac  field  is  assumed  to  be 
uniform.  Substituting  M  =  where  x  is  the  volume  sus¬ 
ceptibility  (dimensionless),  H  =  H0  cos  2ir  ft  is  the  applied 
ac  magnetic  field,  and  rearranging  Eq.  (5),  the  final  result  is 
obtained: 


L'llirfH _  ’ 


(6) 

_l 


where /is  in  hertz  and  x  is  in  SI  units.  To  convert  to  cgs  units, 
the  value  of  ^  is  divided  by  4  tr.  The  assumption  of  uniform  M 
and  H  is  a  good  one  for  relatively  long  cylinders  character¬ 
ized  by  small  demagnetiztion  factors  (see  Sec.  I  C). 

B.  Mutual  Inductance  ot  coaxial  solenoids 

The  mutual  inductance  [Ll2  in  Eq.  (3)]  may  be  calculat¬ 
ed  using  computer  techniques  for  the  general  case  of  coaxial, 
but  not  necessarily  concentric,  thick  solenoids.  Figure  2 
shows  the  geometry  and  defines  some  of  the  variables.  The 
mutual  inductance  is  calculated  by  integrating  A,  due  to  a 
current  density  J  in  coil  1 ,  around  the  turns  of  coil  2.  The 
vector  potential  A  has  only  a  <h  component.  At  the  point 
(r,^,z)  in  cylindrical  coordinates,4 


A*(r,  daaj  dflcos#ln[£  ■+■  (£2  +  r2  +  a2  —  larcos  9 )1/J]|f; 

=  tsL  ("‘da a1  Cdel _ \  = 

4ir  J,„  Jo  \(a2  +  r2  —  larcos  9  )(£2  +  r2  +  a2  —  2arcos  d)'n ){'  4t r 


(7) 


where  £  =  z  —  A.,  f ,  =  z  —  /  /2,  =  2  +  l  /2.  and  the  cylin¬ 

drical  coordinates  of  each  current  element  are  a,  8,  A.  J  is 
assumed  constant  and  is  simply 


J  = 


N1 

/(0.2-a.i) 


(8) 


From  the  expression  forZ.,2  in  Eq.  (3),  and  using  Eqs.  (7)  and 

(8) 


Ln  =  - -  <W,  z)-ds,  (9) 

4tr/(j|2-«n)  J 

where  /  is  canceled.  For  the  purpose  at  hand, 
a l2  2: a, ,  —d/ 2,  where  d  is  the  specimen  diameter.  All  that 
remains  is  to  compute  the  contour  integral  of  A. 

For  purposes  of  computation,  coil  2  is  divided  into 
n  x  m  circular  rings,  usually  not  coincident  with  the  actual 
windings  (see  Fig.  2).  The  radius  of  each  ring  is 

=021  +  («'  -  J)(n2j  -«2i)/n.  ‘  =  1.2 . n,  (10) 

and  its  position  is 

zl=ze~  Hi  +  U  -  J) J  =  l,2,...,m.  (11) 

The  double  integral  A  [given  in  Eq.  |7)]  is  evaluated  for  each 
ring  located  at  (r, ,  zt )  using  a  standard  library  routine.  Then 
the  contour  integral  of  A  in  Eq.  (9)  is  simply 

<W.  zydt  =  —  2  S  2trr,  A{r„z,),  (12) 

J  ntn  i  - 1  j  *  i 

where  dt  is  r  d+.  The  factor  N2/mn  is  introduced  to  scale  L ,  2 
to  reflect  the  actual  number  of  turns  in  coil  2.  If  the  coil 
displacement  zd  is  zero  (the  usual  case  when  the  specimen 
cylinder  is  centered  in  the  pickup  coil),  the  index  j  can  run 
from  1  to  m/2  and  the  result  for  L ,  2  doubled,  thereby  taking 
advantage  of  the  axial  symmetry.  The  number  of  nxm  dis¬ 
crete  segments  should  be  chosen  to  optimize  the  tradeoff 
between  accuracy  and  computation  time. 

Other  procedures  may  be  used  to  calculate  L>2.  Thin- 
coil  approximations  are  clearly  suitable  for  the  model  sole¬ 
noid  and  are  often  adequate  for  the  pickup  coil.5 


I - 

C.  Demagnetization  factor  considerations 

When  susceptibilities  are  large  or  when  cylinders  are 
not  too  long,  the  measured  x  should  be  corrected  for  demag¬ 
netization  fields 

Xm-XA\-DxU  (13) 

where  is  the  "internal"  susceptibility,  corrected  for  de¬ 
magnetization,  and  D  is  the  demagnetization  factor  (0  <  D 
<  I  in  SI).  Exact  values  of  D  are  obtainable  only  for  ellip¬ 
soids.  Average  values  for  uniformly  magnetized  cylinders 
have  been  computed  by  Brownh-7  and  Crabtree"  and  are  giv¬ 
en  in  Table  I.  These  “magnetometric"  values  are  greater 
than  those  "ballistic"  values  reported  by  Bozorth  and  Cha¬ 
pin9  for  highly  permeable  rods.  The  latter  were  experimen¬ 
tally  determined  and  based  on  the  magnetization  at  the  mid¬ 
dle  of  a  rod  rather  than  an  average  over  the  entire  volume. 
Note  that  the  approximation  of  uniform  longitudinal  mag- 
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netization  is  realistic  only  in  the  cases  of  materials  with  small 
susceptibilities,  or  relatively  long  cylinders  characterized  by 
small  demagnetization  factors.  Therefore,  calculations  with 
cylinders  are  always  approximate  to  some  degree. 

D.  Eddy  currant  consideration* 

The  eddy  currents  induced  in  a  conducting  material 
exposed  to  an  alternating  magnetic  field  are  a  function  of  the 
material’s  susceptibility,  resistivity,  and  geometry.  These 
currents  give  rise  to  an  internal  field  which  differs  from  the 
applied  field  in  magnitude  and  phase.  For  accurate  suscepti¬ 
bility  measurements,  these  deviations  should  be  minimized. 

Zijlstra10  examines  the  case  of  an  infinitely  long  cylin¬ 
der  subject  to  an  ac  axial  field.  The  deviations  in  field  magni¬ 
tude  and  phase  are  a  function  of  the  reduced  radius 

a0  =  d/yl2S,  (14) 

where  d  is  the  cylinder  diameter  and  S  is  the  skin  depth, 

6  =  [p/irfp)ul,  (15) 

where ^  is  the  material  resistivity, /is  the  frequency,  and/4  is 
the  permeability  equal  top^X  +  x  I,  where  pQ  is  the  perme¬ 
ability  of  free  space. 

The  analysis  shows  that  magnitude  errors  in  x  are  less 
than  3%  when  a0  <  0.9  and  that  errors  in  both  magnitude 
and  phase  are  negligible  when  a0  <  0.2.  The  restriction  on  a0 
may  be  met  by  reducing  d  and  f,  since  p  and  p  are  materials 
properties. 

II.  EXPERIMENT 
A.  Example 

The  dependence  of  L  *  upon  a  cylinder’s  length  and  di¬ 
ameter  was  computed  for  a  specific  pickup  coil  for  illustra¬ 
tive  purposes.  A  plot  is  shown  in  Fig.  3.  The  normalization  of 
the  variables  applies  to  this  particular  coil  only. 

The  room  temperature  susceptibility  of  a  sample  of 
AlSI-type  316  stainless  steel  was  measured  using  this  coil. 

Table  I.  Longitudinal  demagnetization  factors  D  for  cylinders  as  a  func¬ 
tion  of  the  ratio  of  length  to  diameter  l/d  (after  Crabtree  (Ref.  8)]. 


l/d 

D 

i/d 

D 

I/d 

D 

00 

1000 

2.0 

0181 

4.0 

0.0978 

0.1 

0.796 

2.1 

0.174 

4.1 

0.0956 

0.2 

06*0 

2.2 

0.167 

4.2 

0.0935 

0.3 

0594 

2.3 

0.161 

4.3 

0.0914 

0.4 

0.528 

2.4 

0155 

4.4 

0.0895 

0.3 

0.474 

2.3 

0.149 

4.3 

0.0876 

06 

0.430 

2.6 

0.144 

4.6 

0.0858 

0.7 

0.393 

2.7 

0140 

47 

00841 

OS 

0361 

2.8 

0.135 

4.8 

0.0824 

09 

0334 

2.9 

0.131 

49 

00808 

1.0 

0311 

3.0 

0127 

5.0 

0.0793 

1.1 

0.291 

3  1 

0123 

5  5 

00723 

1.2 

0.273 

3.2 

0120 

6.0 

0.0666 

1.3 

0.257 

3.3 

0116 

6.5 

0.0616 

1.4 

0.242 

34 

0113 

70 

00573 

1.5 

0230 

3  5 

0.110 

7.5 

00536 

16 

0.218 

36 

0107 

80 

00503 

1.7 

0207 

3.7 

0105 

8.5 

00473 

IS 

0.19S 

3.8 

0102 

90 

00447 

19 

0.189 

3.9 

0.100 

100 

0.0403 

The  result  was  verified  with  a  SQUID  susceptometer.  A  cyl¬ 
inder  of  length  10.2  cm  and  diameter  1.2  cm  was  prepared. 
L  •  was  computed  to  be  14/iH  for  the  combination  of  cylin¬ 
der  and  pickup  coil.  An  ac  field  of  570- A/m  rms  at  100  Hz 
was  applied,  x  was  found  to  be  2.9  X  10“ 3  (SI).  This  approxi¬ 
mates  the  smallest  values  of  susceptibility  which  can  be  ac¬ 
curately  measured  with  the  present  apparatus.  A  small  piece 
of  steel  was  obtained  from  the  same  rod  stock  for  SQUID 
measurement.  Using  a  density  of  7.958  g/cm3,"  \  was 
3.0X  10-3,  in  good  agreement  with  the  mutual  inductance 
method.  Taking  a  resistivity  of  75 pfi  cm, 1 2  and  applying  the 
eddy-current  criterion  of  Sec.  I  D,  a0  is  found  to  be  0. 195, 
which  is  very  acceptable. 

B.  Coll  compensation 

As  mentioned  in  Sec.  II  A,  the  background  emf  of  the 
pickup  coil  with  no  sample  present  needs  to  be  subtracted 
from  the  measured  emf  with  the  sample  in  place.  This  may  be 
done  by  the  various  methods  shown  in  Table  II.  Methods  (4), 
(5),  and  (6)  seem  to  be  satisfactory;  method  (6),  illustrated  in 
Fig.  1,  it  the  least  expensive  and  most  straightforward. 

III.  OTHER  CALIBRATION  METHODS 
A.  Chang*  In  Mtf-irviuctanc* 

A  variation  of  the  method  described  in  Sec.  I  A  is  that  of 
relating  the  change  in  self-inductance  of  a  pickup  coil,  after  a 
specimen  is  inserted,  to  the  magnetic  susceptibility  of  the 
specimen.  This  technique  may  be  used  with  a  commercial 
inductance  meter  equipped  with  a  convenient  zero  offset. 
Such  devices  may  be  thought  of  as  applying  a  sinusoidal 
current  to  a  test  coil  and  measuring  the  voltage  across  it.  The 
following  equations  apply: 

(16) 

dt 

where  L,  is  the  self-inductance, 

/  =  /„  cos  hr  ft,  (17) 

t',.™  =2ir  (18) 
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Table  II  Virion  mtthtxh  of  pickup  coil  compenutioa. 


Method  of  coil  compnmrion 


Remarks 


ID 

(2) 

(3) 

(4) 
U) 

(«) 


Subtract  dc  output  of  lock-io  amptifier  for  empty  pickup  coil 

Null  dc  output  of  lock-in  impltfler  using  internal  offset 

Subtract  ac  signal  of  pickup  coil  using  digital  recording  oscilloscope 

Subtract  tins  level  of  ac  signal  of  pickup  coil  using  digital  oscilloscope 

Subtract  an  ac  waveform  equal  to  that  of  empty  pickup  coil,  using  a 

differential  amplifier,  before  input  to  lock-in  amplifier 

Scrim  matched,  counterwound  “bucking”  coil  to  null  sample  pickup  coil 


Simple;  poor  resolution  because  dc/i>  is  small 

Simple;  offset  works  only  over  limited  range 

Any  phase  shift  caused  by  eddy  currents  gives  errors 

Digital  processing  and  noise  filtering  needed 

Need  stable  phase-lock  phase  shifter  with  variable  gain; 

eddy-current  phase  shifts  give  errors 

Both  coils  must  be  exposed  to  same  ac  field 


From  Eq.  (6),  the  sample  contribution  to  the  voltage  is 
—  L*lxhrfHmi.  (19) 

The  field  is  proportional  to  the  current  in  the  coil 

(20) 

where  c  is  a  constant  of  dimension  m~'.  The  total  voltage  is 
0..™,  +  Urm,  =  +  cL  •ly) 

=2i rfl^L.+AL).  (21) 

The  volume  susceptibility  is  obtained  from  the  increase  in 
inductance  AL  when  the  specimen  is  inserted  in  the  coil 


This  method,  relying  on  the  coil  for  both  field  and  pickup,  is 
less  sensitive  than  the  method  in  Sec.  I  A. 

B.  Calibration  with  standards 

It  is  often  desirable  to  calibrate  susceptometers  experi¬ 
mentally,  which  requires  the  use  of  standards.  The  resulting 
calibration  is  strictly  valid  only  for  specimens  of  the  same 
size  and  shape  as  the  standard  used. 

1.  Materials  with  known  susceptibility 

The  NBS  Office  of  Standard  Reference  Materials  sells 
four  susceptibility  standards:  Al,  Pt,  Pd,  and  MnF2,  with  * 
ranging  from  2x  10~5  to  6X  10-3  (SI)  at  room  tempera¬ 
ture.13  They  are  available  in  various  forms. 

2.  Soft  ferromagnets  with  "infinite"  internal  susceptibility 
and  known  demagnetization  factor 

An  example  is  a  spherical  sample  of  Gd  near  its  Curie 
temperature.  Some  provision  for  cooling  is  needed  to  experi¬ 
mentally  obtain  the  maximum  pickup  voltage. 

The  susceptibility  y  is  generally  proportional  to  the 
pickup  coil  voltage  v,  and  inversely  proportional  to  the  sam¬ 
ple  volume  V,  the  ac  field  H,  and  its  frequency /: 

y  =  av/VfH.  (23) 

The  empirical  calibration  constant  a  is  obtained  by  measur¬ 
ing  the  pickup  coil  voltage  t>0  for  a  sphere  of  volume  V0  and 
demagnetization  factor  1/3  in  a  field  H0  of  frequency  f0. 
Assuming  an  infinite  ym ,  Eq.  ( 1 3)  yields 

I  ~/>*o  =  0.  (24) 

whence  the  measured  susceptibility  *0  =  3.  Thus, 

a  =  3F0/0W0/Uo.  (25) 


3.  Superconducting  materials  with  known 
demagnetization  factor 

An  example  is  a  spherical  sample  of  Nb  below  its  criti¬ 
cal  temperature.  The  ability  to  cool  to  cryogenic  tempera¬ 
tures  is  required.  The  perfect  diamagnetism  (*,„,  =  —  1 )  of 
superconductors  arises  from  the  Meissner  effect.  From  Eq. 
(13), 

To  =  DXo  ~  1.  (26) 

and  for  the  sphere,  *0  =  —  3/2,  or  3/2  but  1 80“  out  of  phase. 
Equation  (23)  applies  with 

a  =  3V0f0H0/2v0.  127) 

Similar  calibrations  may  be  done,  of  course,  with  cylin¬ 
drical  standards  using  the  approximate  demagnetization 
factors  in  Table  I.  Methods  of  preparing  spheres  are  de¬ 
scribed  in  Ref.  14. 
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QUENCH  CIRCUIT  FOR  ELECTRONIC  INSTRUMENTS  USED  WITH  SUPERCONDUCTING  MAGNETS 


R.  G.  Benson,  R,  B.  Goldfarb,  and  E.  S.  Pittman 
Electromagnetic  Technology  Division 
National  Bureau  of  Standards 
Boulder,  Colorado  80303 


ABSTRACT 

A  multifunction  circuit  is  described  that  protects,  in  the  event  of  a 
quench,  instruments  connected  or  coupled  to  a  superconducting  magnet. 

INTRODUCTION 

Most  laboratory  superconducting  magnets  are  suitably  protected  against 
damage  if  they  revert  to  the  normal  state  (quench).  Often  it  is  desirable  to 
provide  protection  for  instruments  connected  or  coupled  to  such  a  magnet. 

Such  protection  is  provided  by  a  circuit  that  performs • three  functions:  (1) 
sends  an  interrupt  command  to  a  computer  via  an  interface  card,  (2)  protects 
an  analog-to-digital  (A/D)  card  that  monitors  the  magnet  voltage  from  large 
voltages  generated  during  a  quench,  and  (3)  shorts  the  output  of  a  digital-to- 
analog  (D/A)  programming-voltage  card  and  thus  nulls  the  output  of  the  magnet 
current  supply.  A  diode  bank  that  protects  a  magnet  power  supply  is  also 
described. 

QUENCH  DETECTION  CIRCUIT 

A  superconducting  magnet  in  our  laboratory  is  powered  by  a  constant- 
current  supply.  This  power  supply,  programmed  with  the  dc  voltage  output  of  a 
D/A  card,  is  bipolar  and  referenced  to  ground.  The  D/A  card  is  controlled  by 
a  microcomputer.  The  magnet  voltage,  L  dl/dt,  is  sensed  by  an  A/D  card  which 
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in  turn  is  connected  to  the  computer.  By  monitoring  this  voltage,  the 
computer  is  able  to  determine  when  the  magnet  has  finished  charging  or 
discharging. 


In  the  event  of  a  quench,  we  wish  to:  (1)  inform  the  computer  that  a 
quench  has  occurred  by  enabling  the  interrupt  card,  (2)  protect  the  A/D  card 
against  possible  high  magnet  voltages,  and  (3)  bring  the  power  supply  output 
to  zero  by  nulling  the  programming  voltage  of  the  current  supply. 

We  designed  and  constructed  a  circuit  to  perform  these  tasks.  The  input 
is  connected  across  the  magnet.  The  circuit's  input  impedance  is  greater  than 
500  k 0.  A  comparator-threshold  adjustment  sets  the  voltage  criterion  for  a 
quench  (e.g.,  3.5  V).  A  time-constant  adjustment  prevents  the  quench  detector 
from  triggering  on  random,  transient  voltage  spikes.  This  sets  a  time  delay 
of,  typically,  10  to  250  ms  for  actual  quench  detection.  When  a  quench  is 
detected,  an  LED  turns  on.  Once  the  magnet  voltage  has  decreased  to  below  the 
quench-detection  threshold,  pressing  the  RESET  button  rearms  the  circuit. 

The  circuit  schematic  is  shown  in  Fig.  1.  There  are  three  outputs: 
INTERRUPT,  A/D  OUT,  and  D/A  OUT.  (1)  When  a  quench  is  detected,  a  logical  "1" 
is  sent  to  the  interrupt  card.  (2)  A/D  OUT  normally  sends  the  magnet  voltage 
to  the  A/D  card.  In  the  event  of  a  quench,  a  relay  breaks  this  connection. 
Additional  protection  is  provided  by  a  series  fuse.  Finally,  two  Zener 
diodes,  D3,  in  parallel  with  the  output  prevent  the  magnitude  of  the  output 
voltage  from  exceeding  12  V.  (3)  D/A  OUT  is  normally  an  open  switch  in 
parallel  with  the  output  of  the  D/A  card.  When  a  quench  occurs.  Relay  A 
closes,  forcing  the  programming  voltage  of  the  magnet  current  supply  to  zero. 
(The  D/A  card  is  internally  current  limited.) 

The  circuit  is  divided  into  several  stages:  input  protect,  absolute 
value,  voltage  comparator,  latch  logic,  and  output  stages. 
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Input  protect.  Large  voltage  spikes  dissipate  power  in  resistor  Rl.  The 
input  voltage  to  the  rest  of  the  circuit  is  limited  to  about  10  V  by  the  two 


Zener  diodes  Dl. 

Absolute  value.  The  absolute  value  stage  uses  two  operational  ampli¬ 
fiers,  A1  and  A2,  and  has  a  positive  output  with  unity  gain.  For  positive 
input  voltages,  diode  D2  is  reverse  biased,  A2  is  effectively  removed  from 
the  circuit,  and  A1  is  a  unity-gain  amplifier.  For  negative  voltages,  D2  is 
forward  biased,  the  positive  input  to  M  is  grounded,  and  A1  acts  as  an 
inverter. 

Voltage  comparator.  The  comparator  produces  a  low  output  of  zero  volts 
for  inputs  greater  than  the  reference  voltage  and  a  high  output  of  4  V  other¬ 
wise.  The  feedback  resistor  provides  some  hysteresis  in  the  comparator.  The 
QUENCH-THRESHOLD  adjustment  R2  sets  the  comparator  reference  voltage  V 
The  TIME-CONSTANT  adjustment  R3  is  part  of  a  simple  RC  filter. 

Latch  logic.  When  the  circuit  is  turned  on,  the  10-pF  capacitor  Cl  holds 
the  CLEAR  input  to  ground  long  enough  to  allow  the  flip-flop  to  start  in  the 
cleared  state  (Q  low  and  Q  high).  After  about  10  ms,  the  CLEAR  input  goes 
high.  It  can  be  set  low  with  the  RESET  button.  When  a  quench  is  detected, 
the  output  of  the  voltage  comparator  goes  low  and  SET  becomes  active,  causing 
Q  to  go  high  and  Q  low. 

Output  stages.  When  Q  is  high,  a  logical  "1"  is  sent  to  the  interrupt 
card.  Also,  transistor  Tl  is  turned  on,  which  lights  the  LED  and  closes  Relay 

A.  This  in  turn  shorts  D/A  OUT,  which  zeros  the  programming  voltage  to  the 
power  supply.  When  Q  is  low,  transistor  T2  is  turned  off.  This  opens  Relay 

B,  which  protects  the  A/D  card. 
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POWER  SUPPLY  PROTECTION 


Power  supplies  designed  for  use  with  superconducting  magnets  are  often 
built  with  quench  protection.  A  quench  is  accompanied  by  large  voltages 
(typically  a  few  hundred  volts)  and  currents  of  polarity  opposite  that  of  the 
charge  current.  Most  power  supplies  are  not  designed  to  sink  reverse  current. 
In  the  case  of  a  unipolar  dc  supply,  adequate  protection  may  be  provided  by  a 
reverse-biased,  heat-sinked  power  diode  across  the  output  terminals,  in 
parallel  with  the  magnet  leads. 

In  the  case  of  a  bipolar  supply,  however,  the  forward  voltage  drop  of  the 
diode  would  be  too  small  for  negative-current  operation;  magnet  charging 
voltages  are  often  as  great  as  3  V.  We  have  constructed  an  array  of  heat- 
sinked  1N3111  diodes  to  protect  such  bipolar  supplies.  The  circuit  is  shown 
in  Fig.  2.  It  is  connected  across  the  output  of  the  power  supply,  in  parallel 
with  the  magnet  leads.  The  diode  array  does  not  conduct  for  typical  magnet 
charging  voltages.  If  quench  voltages  are  large  in  either  polarity,  the  array 
conducts,  thereby  protecting  the  power  supply.  The  six  diodes  in  series  in 
either  direction  begin  to  conduct  at  4  V  and  conduct  well  at  5  V. 
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I.  PHYSICAL  CHARACTERISTICS  OF  MULTIFILAMENTARY  Nb-Ti  SUPERCONDUCTORS 
We  have  found  that,  except  for  the  number  of  Nb-Ti  filaments,  super¬ 
conductor  wire  manufacturers  often  do  not  reliably  report  the  physical  char¬ 
acteristics  of  their  wires.  For  accurate  magnetization  and  loss  calculations, 
however,  we  must  know  the  copper-to-superconductor  volume  ratio  and  the 
filament  radius  in  addition  to  measured  magnetic  moment.  It  is  also  useful  to 
know  the  density  of  the  Nb-Ti  alloy.  This  note  describes  how  we  determine 
these  parameters. 

A  1.000-meter  length  of  wire  is  stripped  of  its  insulation.  We  measure 
its  mass,  mwire-  Using  Archimedes'  principle,  we  determine  its  volume,  vwire- 
The  wire  is  then  soaked  in  a  50%  solution  of  nitric  acid  and  water  to  remove 
the  copper,  leaving  only  the  Nb-Ti  filaments.  The  filaments  are  rinsed  and 
thoroughly  dried.  We  then  measure  their  mass,  m^^.  We  do  not  directly 
determine  the  volume  of  the  filaments,  because  of  the  danger  of  trapping 

air  bubbles  upon  immersion  in  water.  Instead,  we  calculate  the  mass  of 
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copper,  mCu  ■=  mwire  “  mfii*  anc*  use  t^ie  ^nown  density  of  pure  Cu  (8.93  g/cin  ) 

to  calculate  the  volume  of  copper,  V.  .  Then,  V,..  *  V  .  -  V_  . 

Cu  fil  wire  Cu 

The  copper-to-superconductor  volume  ratio  is  then  V^/V^^.  The  density 

U 

of  the  Nb-Ti  alloy  is  The  filament  radius  is  (V^^/irNL)  ,  where  N 

is  the  number  of  filaments  given  by  the  manufacturer  and  L  is  the  length 

(1.000  m). 
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II.  CONVERSION  OF  ATOMIC  AND  WEIGHT  PERCENTS 
Let  the  atomic  masses  of  the  elements  be  denoted  by  A^,  the  atomic 
percents  be  denoted  by  a^,  and  the  weight  percents  be  denoted  by  w^.  By 
definition,  Ea^  =  1  and  Ew^  =  1. 

To  convert  atomic  to  weight  percent: 

Wi  *  aiAi/EajAj  •  (1) 

The  average  atomic  mass  of  the  alloy  is  A: 

A  =  Ea  A  .  (2) 

Thus, 

wA  *  a^A^/A  .  (3) 

To  convert  weight  to  atomic  percent: 

ai  ”  (wi / Ai> /E (w . /A . )  .  (4) 

As  pointed  out  by  H.  M.  Ledbetter  (private  communication),  substituting  Eq.  4 
into  ^q.  2  gives 

A  =  [E(w  /A  n”1  .  (5) 

Thus, 

ai  *  (w1/Ai)  A  .  (6) 

III.  DENSITIES  OF  SUBSTITUTIONAL  BINARY  ALLOYS 
Often  we  wish  to  deduce  the  density  of  a  binary  alloy  by  interpolation 
from  the  densities  of  the  pure  elements.  There  are  two  approaches.  The  first 
is  to  assume  that  both  elements  have  the  same  atomic  volumes.  We  then 
interpolate  by  substituting  the  atomic  mass  of  one  for  the  atomic  mass  of  the 
other.  The  second  approach  is  to  assume  the  same  atomic  masses.  We  then 
interpolate  by  substituting  the  atomic  volumes.  Both  methods  are  presented 
below,  applied  to  the  Nb-Ti  system.  In  this  example,  where  the  atomic  masses 
of  Nb  and  Ti  are  known  to  be  very  different,  the  first  approach  is  better. 


Subscripts  indicate  element  1  or  2 
p  *  density 
m  *  mass  per  atom 
v  *  volume  per  atom 
n  *  number  of  atoms  per  cm3 

Same  Atomic  Volumes 


q  =  number  of  atoms  per  g 
a  =  atomic  percent  element  1 
w  ■  weight  percent  element  1 


There  is  the  same  number  of  atoms  per  cm3  for  elements  1  and  2  (n^  *  n^) . 
The  atomic  mass  of  one  is  substituted  for  the  atomic  mass  of  the  other. 

p  =  nm  (7) 

Palloy  *  n[aml  +  <1"a)m21 

*  ap1  +  (l-a)p2  (8) 

For  example,  PNfa  *  8.57  g/cm3  and  p,^  »  4.51  g/cm3.  A  "standard"  alloy 

composition  is  46.5  wt.%  Ti  *  63.0  at.%  Ti.  Thus, 

PNbTi  **  6,012  8/cm3  (9) 

This  compares  favorably  with  a  measured  density  of  6.004  g/cm3. 

Same  Atomic  Masses 

There  is  the  same  number  of  atoms  per  g  for  elements  1  and  2  (q^  =  q2> . 
The  atomic  volume  of  one  is  substituted  for  the  atomic  volume  of  the  other. 


Palloy  *  9[UV1  +  <1"u>v21 

*  u/p1  +  ( 1— u) / p £  (11) 

Palloy  *  PiP2/[(pl+u(p2"Pl)]  (12) 

For  the  example  of  Nb-Ti  of  "standard"  composition,  we  would  obtain  an  inac¬ 


curate  estimate:  p. 


5.784  g/cm3.  The  first  method,  assuming  equal 


atomic  volumes  for  the  elements,  is  more  suitable. 

If  the  densities  of  certain  alloy  compositions  are  known,  further  refine¬ 
ment  is  possible  by  interpolating  between  them  instead  of  the  densities  of  the 
pure  elements. 
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UNITS  FOR  MAGNETIC  PROPERTIES 


Quantity 

Symbol 

Gaussian  &  cgs  emu' 

Conversion 
factor.  C  b 

SI  &  rationalized  mks  ‘ 

Magnetic  Dux  density, 
magnetic  induction 

B 

gauss  (G)  d 

10  4 

tesla  (T),  Wb/m; 

Magnetic  flux 

d> 

maxwell  (Mx).  G-cnr 

10  * 

weber  (Wb).  volt  second  (V  s) 

Magnetic  potential  difference, 
magnetomotive  force 

U.  F 

gilbert  (Gb) 

10/47T 

ampere  (A) 

Magnetic  field  strength, 
magnetizing  force 

H 

oersted  (Oe)/  Gb/cm 

IO'/4ir 

A/m  ^ 

(Volume)  magnetization  8 

U 

emu/cm!  * 

10' 

A/m 

(Volume)  magnetization 

4  irM 

G 

10’/4tr 

A/m 

Magnetic  polarization, 
intensity  of  magnetization 

J.  1 

emu/cm1 

irr  X  10-4 

T.  Wb/m! ' 

(Mass)  magnetization 

<T.  M 

emu/g 

1 

4rr  X  10 

Anr/kg 

Wbm/kg 

Magnetic  moment 

m 

emu,  erg/G 

10"' 

Anr.  joule  per  tesla  (J/T) 

Magnetic  dipole  moment 

j 

emu,  erg/G 

4irx  10- 10 

Wb  m ' 

(Volume)  susceptibility 

X ■  * 

dimensionless,  emu/cm1 

4tt 

(4#)^  10-7 

dimensionless 

henry  per  meter  (H/m),  Wb/(A  m) 

(Mass)  susceptibility 

Xe •  *e 

cmVg,  emu/g 

4rrx  10-’ 
(4ir):x  10' 10 

mJ/kg 

Hm;/kg 

(Molar)  susceptibility 

^molr  tfmol 

cmVmol,  emu/mol 

4?rx  lO'1’ 
(47r)Jx  10" 11 

mVmol 

H-m’/mol 

Permeability 

d 

dimensionless 

4  jx  10'’ 

H/m.  Wb/(A  m) 

Relative  permeability7 

Pr 

not  defined 

dimensionless 

(Volume)  energy  density, 
energy  product  * 

w 

erg/cm1 

10'1 

J/mJ 

Demagnetization  factor 

D,  .V 

dimensionless 

l/4rr 

dimensionless 

a.  Gaussian  units  and  cgs  emu  are  the  same  for  magnetic  properties.  The  defining  relation  is  B  =H  +  iw\f . 

b.  Multiply  a  number  in  Gaussian  units  by  C  to  convert  it  to  SI  (e.g„  1  Gx  10-4  T/G  =  10"4  T). 

c.  SI  (Systeme  International  d'Unites)  has  been  adopted  by  the  National  Bureau  of  Standards.  Where  two  conversion  factors  are 
given,  the  upper  one  is  recognized  under,  or  consistent  with.  SI  and  is  based  on  the  definition  B  =  po(ff  ■*-  M ),  where 
fio  =  4rrx  10'"  H/m.  The  lower  one  is  not  recognized  under  SI  and  is  based  on  the  definition  B  =  naff  W,  where  the  symbol 
/  is  often  used  in  place  of  J 

d.  1  gauss  =  10!gamma(y). 

e.  Both  oersted  and  gauss  are  expressed  as  cm  *,/J  gl/J-s''  in  terms  of  base  units. 

/  A/m  was  often  expressed  as  “ampere-turn  per  meter”  when  used  for  magnetic  field  strength. 

g.  Magnetic  moment  per  unit  volume, 

h.  The  designation  "emu"  is  not  a  unit. 

/.  Recognized  under  SI.  even  though  based  on  the  definition  B  =  po/f  +J.  See  footnote  c. 

j.  pr  =  p/po  —  1  +X-  all  in  SI-  dr  is  equal  to  Gaussian  p. 

k.  BH  and  po M  H  have  SI  units  J/m1;  M-H  and  BH/ 4rr  have  Gaussian  units  erg/cm1. 

R.  B.  Goldfarb  and  F  R.  Fickett.  U  S.  Department  of  Commerce,  National  Bureau  of  Standards.  Boulder.  Colorado  80303,  March  1985 
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